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Pulse  Propagation  and  Bistatic  Scattering 


Abstract 

A  set  of  pulse-propagation  coupling  equations  is  derived.  They  couple  the  output  electrical 
signal  at  a  point  element  in  a  receive  array  to  the  transmitted  electrical  signal  at  the  input  to  a 
transmit  array  via  the  complex  frequency  response  of  a  fluid  medium  (e.g.,  air  or  water).  The  pulse- 
propagation  coupling  equations  are  based  on  linear,  time-variant,  space-variant,  filter  theory,  the 
principles  of  complex  aperture  theory  and  array  theory,  and  solving  a  linear  wave  equation,  which 
includes  satisfying  all  boundary  conditions,  including  the  boundary  condition  at  the  source.  They 
can  be  used  to  accurately  model  the  propagation  of  small-amplitude  acoustic  pulses  in  the  ocean  for 
a  bistatic  scattering  problem.  Three  different  bistatic  scattering  problems  are  considered:  1)  no 
motion,  2)  only  the  discrete  point  scatterer  is  in  motion,  and  3)  all  three  platforms  (the  transmitter, 
discrete  point  scatterer,  and  receiver)  are  in  motion.  Specific  examples  on  the  use  of  the  pulse- 
propagation  coupling  equations  are  given  for  the  three  different  bistatic  scattering  problems. 


1  Introduction 

The  main  purpose  of  this  report  is  to  establish  a  mathematical  framework  that  can  be  used  to 
accurately  model  the  propagation  of  small-amplitude  acoustic  pulses  in  the  ocean  for  a  bistatic 
scattering  problem.  Since  we  are  dealing  with  small-amplitude  acoustic  pulses,  a  linear  wave 
equation  accurately  describes  the  propagation  of  sound  between  source,  discrete  point  scatterer,  and 
receiver.  As  a  result,  we  treat  the  propagation  of  small-amplitude  acoustic  pulses  in  the  ocean  as 
transmission  through  a  linear,  time-variant,  space-variant  filter.  Treating  the  ocean  medium  as  a 
linear  filter  is  valid  because  we  are  trying  to  solve  a  linear  wave  equation. 

Part  2  of  this  report  provides  necessary  background  information  on  some  of  the  principles  of 
linear,  time-variant,  space-variant  filter  theory  and  how  it  relates  to  solving  a  linear  wave  equation.  In 
particular.  Section  2.1  introduces  the  time-variant,  space-variant,  complex  frequency  response  of 
the  ocean,  which  is  shown  to  be  the  time-harmonic  solution  of  a  linear  wave  equation  when  the 
source  distribution  is  a  time-harmonic  point  source.  Three  different  examples  are  worked  out  in 
Section  2.1. 

Section  2.2  deals  with  the  main  problem  of  pulse  propagation  from  a  linear  systems  theory 
point  of  view.  In  this  section  we  derive  a  set  of  equations  that  we  refer  to  as  the  pulse-propagation 
coupling  equations.  They  couple  the  output  electrical  signal  at  a  point  element  in  a  receive  array  to 
the  transmitted  electrical  signal  at  the  input  to  a  transmit  array  via  the  complex  frequency  response 
of  a  fluid  medium  (e.g.,  air  or  water).  The  pulse-propagation  coupling  equations  are  based  on 
linear,  time-variant,  space-variant,  filter  theory,  the  principles  of  complex  aperture  theory  and 
array  theory,  and  solving  a  linear  wave  equation,  which  includes  satisfying  all  boundary 
conditions,  including  the  boundary  condition  at  the  source.  They  provide  a  consistent,  logical,  and 
straightforward  mathematical  framework  for  the  solution  of  small-amplitude,  acoustic,  pulse- 
propagation  problems.  The  main  features  of  the  pulse-propagation  coupling  equations  are  as 
follows:  1)  transmitted  electrical  signals  are  modeled  as  amplitude-and-angle-modulated  carriers, 
2)  both  the  transmit  and  receive  apertures  are  modeled  as  volume,  conformal  arrays  composed  of 
unevenly-spaced,  complex-weighted,  point  elements  (this  type  of  model  for  both  of  the  apertures 
allows  for  maximum  flexibility),  3)  the  complex  weights  are  frequency  dependent  and  allow  for 
beamforming,  4)  the  performance  of  the  point  elements  in  both  the  transmit  and  receive  arrays  are 
characterized  by  frequency-dependent,  transmitter  and  receiver  sensitivity  functions,  and  5)  the 
solution  of  a  linear  wave  equation  is  given  by  the  complex  frequency  response  of  the  fluid  medium. 
It  is  important  to  note  that  attention  to  all  proper  units  of  measurement  are  taken  into  account  in 
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order  to  ensure  the  accurate  prediction  of  signal  strength  levels  at  each  element  in  a  receive  array. 
This  is  especially  important,  for  example,  in  order  to  obtain  accurate  probability  of  detection  results. 

Part  3  of  this  report  is  devoted  to  the  derivation  of  the  complex  frequency  response  of  the 
ocean  for  three  different  bistatic  scattering  problems.  Scatter  from  a  discrete  point  scatterer  is 
modeled  via  the  scattering  amplitude  function,  which  is  a  complex  function  (magnitude  and  phase) 
and  is,  in  general,  a  function  of  frequency,  the  direction  of  wave  propagation  from  the  source  to  the 
scatterer,  and  the  direction  of  wave  propagation  from  the  scatterer  to  the  receiver.  In  addition  to  the 
scattering  amplitude  function,  frequency-dependent  attenuation  is  taken  into  account  in  order  to 
model  the  propagation  of  sound  from  transmitter  to  discrete  point  scatterer,  and  from  discrete  point 
scatterer  to  receiver.  Section  3.1  is  devoted  to  the  first  bistatic  scattering  problem,  which  involves  no 
motion  -  the  transmitter,  the  discrete  point  scatterer,  and  the  receiver  are  not  in  motion.  Section  3.2  is 
devoted  to  the  second  bistatic  scattering  problem  where  only  the  discrete  point  scatterer  is  in  motion, 
and  Section  3.3  is  devoted  to  the  third  bistatic  scattering  problem  where  all  three  platforms  are  in 
motion.  The  dimensionless,  time-compression  /  time-stretch  factor  is  discussed  in  Sections  3.2  and 
3.3  due  to  the  consideration  of  motion.  The  time-compression  /  time-stretch  factor  takes  into 
account  the  relativistic  effects  of  motion  and  provides  for  accurate  time  delay  and  Doppler  shift 
values.  In  addition  to  other  examples,  specific  examples  on  the  use  of  the  pulse-propagation 
coupling  equations  are  given  in  all  three  sections  of  Part  3  for  the  three  different  bistatic  scattering 
problems. 

2  The  Ocean  Medium  As  A  Linear,  Time- Variant,  Space- Variant  Filter 
2.1  Complex  Frequency  Response  Of  The  Ocean 

The  propagation  of  small-amplitude  acoustic  signals  in  the  ocean  medium  can  be  described 
by  the  following  linear ,  three-dimensional,  inhomogeneous  wave  equation: 

V2yMfo  r)  — ITT  for)  =  xM(t,r),  t2-1'1) 

c  (r)  at 

where  yM(t,r)  is  the  velocity  potential  at  time  t  and  position  r  =  (x,y,z)  with  units  of  squared- 
meters  per  second,  xM(t,  r)  is  the  source  distribution  at  time  t  and  position  r  with  units  of  inverse 
seconds,  and  c(r)  is  the  speed  of  sound  in  meters  per  second,  shown  as  a  function  of  position  r. 
The  source  distribution  represents  the  volume  flow  rate  (source  strength)  per  unit  volume  of  fluid. 
The  solution  of  (2.1-1)  can  be  obtained  by  treating  the  ocean  medium  as  a  linear ,  time-variant, 
space-variant  filter.  Treating  the  ocean  medium  as  a  linear  filter  is  valid  because  we  are  trying  to 
solve  a  linear  wave  equation.  The  time-variant  property  of  the  filter  allows  for  motion  and  the 
possibility  that  the  properties  of  the  ocean  medium  may  change  with  time.  The  space-variant 
property  of  the  filter  allows  for  the  presence  of  boundaries,  both  the  speed  of  sound  and  the 
ambient  (equilibrium)  density  of  the  ocean  medium  as  functions  of  position  (spatial  coordinates), 
and  discrete  point  scatterers. 

Figure  2.1-1  is  an  illustration  of  a  linear,  time-variant,  space-variant,  ocean-medium  filter  with 
input-output  relationship  given  by  [1] 

(2.1-2) 

where  the  source  distribution  xM(t,r)  is  the  input  acoustic  signal  to  the  filter,  ftM(t,r|r0,r0)  is  the 
time-variant,  space-variant,  impulse  response  ( Green’s  function )  of  the  filter,  that  is,  it  is  the 


Figure  2.1-1.  Illustration  of  a  linear,  time-variant,  space-variant,  ocean-medium  filter. 
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response  of  the  filter  at  time  t  and  position  r  =  (x,y,z)  due  to  the  application  of  a  unit-amplitude 
impulse  at  time  t0  and  position  r0  =  (x0,y0,z0),  the  velocity  potential  yM(t, r)  is  the  output  acoustic 
signal  from  the  filter,  and  dr0  =  dx0dy0dz0.  Therefore,  the  right-hand  side  of  (2.1-2)  is  shorthand 
notation  for  a  four-fold  integral.  Note  that  (2.1-2)  is  not  a  multidimensional  convolution  integral. 
Equation  (2.1-2)  is  the  solution  of  the  linear  wave  equation  given  by  (2.1-1). 

Also  shown  in  Fig.  2.1-1  is  the  time-variant,  space-variant,  transfer  function  HM(t,r\f,  v)  of 
the  ocean-medium  filter  defined  by  [1] 


HM(t,r\f, v)  4  J_~  J~  hM(t,r\t0,r0) cxp[~j2nf{t - t0)]exp[+j2nv(r - r0)]dt0 dr0 
=  j°°J^hM(t,r\t0,r0)exp[-j2nf(t  -  t0)]dt0  exp[+j2nv  (r  -  r0)]^r0 
=  J_~  Hm (t, r  | /, r0 )  txp[+j2nv  •  (r  -  r0 )]  dr0 , 


(2.1-3) 

where  /  represents  input  (transmitted)  frequency  components  in  hertz,  v  =  (/x,/y,/z)  is  a  three- 
dimensional  vector  whose  components  are  input  (transmitted)  spatial  frequencies  in  cycles  per 
meter  in  the  X,  Y,  and  Z  directions,  respectively,  dr0  =  dx0  dy0  dz0,  and 


HM(t,r\f,r0)  =  j^hM(t,r\t0,r0)exp[-j2nf(t  -  t0)]dt0 


(2.1-4) 


is  the  time-variant,  space-variant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz. 
Similarly  [1], 

hM(t,  r  |  r0,r0)  =  O-L  HM  (t,  r  |  /,  v)  zxp[+j2nf(t  -  t0)]exp[-j2nv  (r  -  r0)]dfdv 
=  J”^J^//A/(t,r|/,v)exp[-;2^v*(r  -  r0)]dvexp[+j2nf(t  -  t0)\df 
=  H M  (t,  r  |  /,  r0 )  exp[+j27r/(?  - 10 )]  df , 

(2.1-5) 


where  dv  =  dfx  dfY  dfz .  From  the  right-hand  side  of  (2.1-5),  it  can  be  seen  that  the  complex 
frequency  response  H M (t, r | f. r0 can  be  obtained  from  the  transfer  function  HM(t,r\f ,vj  as 
follows: 


HM(t,r  |/,r0)  =  j~"HM(t,r\f,  v)exp[-j2nv  (r  -  r0)]rfv. 


(2.1-6) 
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Equation  (2.1-6)  is  the  inverse  of  the  last  relationship  shown  in  (2.1-3). 

The  strategy  that  we  are  going  to  follow  is  to  first  find  time-harmonic  solutions  of  the  linear 
wave  equation  for  various  bistatic  scattering  problems,  and  then,  by  using  Fourier  transform 
techniques,  we  will  obtain  the  corresponding  pulse  solutions.  Recall  that  a  time-harmonic  acoustic 
field  is  one  whose  value  at  any  time  and  point  in  space  depends  on  a  single  frequency  component. 
For  convenience,  let  us  rewrite  (2.1-2)  as  follows: 

=  J  (2-1-7) 

Let  the  source  distribution  xM(t,r)  be  motionless  and  time-harmonic,  that  is, 

xM(t,r)  =  xfM  (r)exp(+j27tft) ,  (2.1-8) 

where  xf  M(r)  is  the  spatial-dependent  part  (with  units  of  inverse  seconds)  of  the  motionless,  time- 
harmonic  source  distribution.  Substituting  (2.1-8)  into  (2.1-7)  yields 

yM(t, r)  =  j  mxf >w(Oj ^hM(t,r\t0£)exp(+j2nft0)dt0d£,  (2.1-9) 

and  upon  multiplying  the  right-hand  side  of  (2.1-9)  by 

exp(-y'2/r/0  exp(+j2nft)  =  1 , 

we  obtain  the  following  expression  for  the  output  acoustic  signal  (velocity  potential)  from  a  linear, 
time-variant,  space-variant,  ocean-medium  filter: 


yM  (t,  r)  =  I _  Xf  -  *  (C)  Hm  (/,  r  I  /,  C)  d&xp(+j2jift) , 


(2.1-10) 


where  |/,r0)  is  the  time-variant,  space-variant,  complex  frequency  response  of  the  ocean  at 

frequency  /  hertz  given  by  (2.1-4). 

If  the  source  distribution  xM(t,  r)  is  a  motionless,  time-harmonic,  point  source  with  units  of 
inverse  seconds,  then  (2.1-8)  is  still  applicable  but 

*/iM(r)  =  S0<5(r-r0),  (2.1-11) 

where  S0  is  the  source  strength  in  cubic  meters  per  second  and  the  impulse  function  <5(r-r0),  with 
units  of  inverse  cubic  meters,  represents  a  point  source  at  r0  =  (x0,y0,z0).  Substituting  (2.1-1 1)  into 
(2.1-10)  yields 


yM(r,r)  =  SQHM(t,r\f,r0)exp(+j27tft), 


(2.1-12) 


where  HM(t,r\f,r0)  is  given  by  (2.1-4). 

Different  forms  of  the  input-output  relationship  given  by  (2.1-2)  and  the  time- variant,  space- 
variant,  complex  frequency  response  given  by  (2.1-4)  can  be  obtained  by  considering  the  following 
three  examples. 
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Example  2.1-1  Time-Invariant,  Space-Invariant,  Ocean-Medium  Filters 

If  the  ocean-medium  filter  is  time-invariant  and  space-invariant,  then  [1] 

hM(t,r\t0,r0)  =  hM(t-t0,r-r0),  (2.1-13) 

that  is,  the  impulse  response  depends  only  on  the  time  difference  t  =  t-t0  and  the  vector  spatial 
difference  R  =  r-  r0.  The  time  difference  T  is  a  measure  of  how  long  ago  the  impulse  was  applied, 
or  the  travel  time  between  source  and  field  points.  The  vector  spatial  difference  R  corresponds  to 
the  distance  and  direction  between  source  and  field  points.  The  magnitude  of  the  vector  spatial 

difference  |R|  =  |r-r0|  is  a  measure  of  how  far  away  the  impulse  was  applied,  or  the  distance 
traveled  between  source  and  field  points.  Substituting  (2.1-13)  into  (2.1-2)  yields 

=  J_m  ~  ^o’** —  ro  )dto  dro  *  (2.1  - 14) 

which  is  a  multidimensional  convolution  integral  as  expected.  In  addition,  substituting  (2.1-13)  into 
(2.1-4)  yields 

Hm (t,r |/,r0)  =  £  hM(t  -t0, r - r0)exp[~72^/(r  - t0)]dt0 .  (2.1-15) 

If  we  let  T  =  t  -  f0.  then  dt0  - -dr  and  (2.1-15)  reduces  to 

HM(t,r\f,r0)  =  £  hM(t,r  —  r0)  dr 

=  FT{hM(r,r-r0)}  (2.1-16) 

=  HM(f,  r-r0). 

If  the  source  distnbution  is  motionless  and  time-harmonic  [see  (2.1-8)]  and  the  ocean- 
medium  filter  is  time -invariant  and  space-invariant,  then  substituting  (2.1-16)  into  (2.1-10)  yields 
the  following  expression  for  the  output  acoustic  signal  (velocity  potential)  from  the  ocean-medium 
filter: 

yM(t,r)  =  y/M(r)exp(+j2nft),  (2.1-17) 

where 

y/JW(r)  =  xfM(r)*HM(f,r)  =  Qd^  (2.1-18) 
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and  the  time-invariant,  space-invariant,  complex  frequency  response  of  the  ocean  at  frequency  / 
hertz  is  given  by  [see  (2.1-16)] 


H„(f,  r  -  r0)  =  Fz{hM(t,  r  -  r0)}  =  \_hM(r,r  -  r0)exp  {-j2nfx)dx. 


(2.1-19) 


Furthermore,  if  the  source  distribution  is  a  motionless,  time-harmonic,  point  source,  then  (2.1-8) 
and  (2.1-11)  are  applicable,  and  upon  substituting  (2.1-11)  into  (2.1-18),  we  obtain 


(2.1-20) 


where  HM(f, r  -  r0)  is  given  by  (2.1-19). 

For  example,  the  impulse  response  of  an  ideal  ( nonviscous ),  unbounded,  homogeneous, 
ocean  medium  is  given  by  [2, 3] 


Mr’rl'o’ro)  hM(t  t0,T  r0) 


4tt  |  r  —  ] 


t-tn 


(2.1-21) 


Equation  (2.1-21)  indicates  that  the  ocean-medium  filter  is  time-invariant  (no  motion  and  the 
properties  of  the  medium  between  source  and  field  points  do  not  change  with  time)  and  space- 
invariant  (unbounded,  homogeneous  medium  with  no  discrete  point  scatterers  between  source  and 
field  points).  If  we  let  r  =  t-t0,  then  substituting  (2.1-21)  into  (2.1-19)  yields  the  following 
expression  for  the  time-invariant,  space-invariant,  complex  frequency  response  of  the  ocean  at 
frequency  /  hertz: 


where 


r-r0) 


exp(-yA:|r-r0|) 

4^|r-r0| 


(2.1-22) 


k  =  Infjc  =  In /X 


(2.1-23) 


is  the  wavenumber  with  units  of  radians  per  meter,  c  is  the  constant  speed  of  sound  of  the 
homogeneous  ocean  medium  in  meters  per  second,  and  X  is  the  wavelength  in  meters.  Note  that 
HM(f, r-r0)  given  by  (2.1-22)  is  an  exact  solution  of  the  following  lossless,  inhomogeneous, 
Helmholtz  equation  [4, 5]: 


V2HM(f,r-r0)  +  k2HM(f,r-r0)  =  d(r-r0).  (2.1-24) 

Also  note  that  hM(t-t0, r-r0)  given  by  (2.1-21)  can  be  used  in  (2.1-14)  and  that  HM(f, r-r0) 
given  by  (2.1-22)  can  be  used  in  both  (2.1-18)  and  (2.1-20),  and  that 
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(2.1-25) 


hm(M  =  - 


exp(-yl|r|) 

An\r\ 


Example  2.1-2  Time-Invariant,  Space- Variant,  Ocean-Medium  Filters 

If  the  ocean-medium  filter  is  time-invariant  and  space-variant ,  then  [1] 


lif^  t0,rQ)-hM(t  ^r|r0), 


(2.1-26) 


that  is,  the  impulse  response  depends  on  the  time  difference  r  =  t  —  t0.  The  time  difference  r  is  a 
measure  of  how  long  ago  the  impulse  was  applied,  or  the  travel  time  between  source  and  field 
points.  Substituting  (2.1-26)  into  (2.1-2)  yields 


y>A (*> r)  —  J  x  J  ^ (t0’ro )K if  Vlro )dto dr0  • 


(2.1-27) 


In  addition,  substituting  (2.1-26)  into  (2.1-4)  yields 


HM(t,r |/,r0)  =  - 10 , r  |  r0 ) exp[-j2nf  (t  -  t0)]dt0. 


(2.1-28) 


If  we  let  r  =  t  - t0,  then  dt0  =  -dr  and  (2.1-28)  reduces  to 


HM(t,r  |/,r0)  =  J_~ K  (t,  r  |  r0 )  exp(-j2nfr)  dr 

=  FT{hM(r,r\r0)} 

=  HM(f,  r|r0). 


(2.1-29) 


If  the  source  distribution  is  motionless  and  time-harmonic  [see  (2.1-8)]  and  the  ocean- 
medium  filter  is  time-invariant  and  space-variant,  then  substituting  (2.1-29)  into  (2.1-10)  yields  the 
following  expression  for  the  output  acoustic  signal  (velocity  potential)  from  the  ocean-medium 
filter: 


yM(t,  r)  =  [jf'M(QHM(  /,  r  I  tyd&xp(+j2nft). 


(2.1-30) 


where  the  time-invariant,  space-variant,  complex  frequency  response  of  the  ocean  at  frequency  / 
hertz  is  given  by  [see(2.1-29)] 


^(/,r|r0)  =  FT{fc„(T,r|r0)}  =  \~~hM{r,r\r0)exV(-j2nfr)dr. 


(2.1-31) 
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Furthermore,  if  the  source  distribution  is  a  motionless ,  time-harmonic,  point  source,  then  (2.1-8) 
and  (2.1-11)  are  applicable,  and  upon  substituting  either  (2.1-11)  into  (2.1-30),  or  (2.1-29)  into 
(2.1-12),  we  obtain 


yj#(*.r)  =  SoHm(  f,r\r0)exp(+j2nft). 


(2.1-32) 


where  HM  ( /,  r  |  r0 )  is  given  by  (2. 1  -3 1 ).  Several  different  examples  of  HM  ( /,  r  |  r0 )  are  given  next. 

A  simple  waveguide  model  of  the  ocean  is  illustrated  in  Fig.  2.1-2  where  p,  c,,  p2c2 ,  and  p3c3 
are  the  characteristic  impedances  of  fluid  media  I  (air),  II  (seawater),  and  III  (ocean  bottom), 
respectively.  The  parameters  p,  and  c;,  i  =  1,2,3,  are  the  constant  ambient  (equilibrium)  densities 
(in  kilograms  per  cubic  meter)  and  speeds  of  sound  (in  meters  per  second)  in  each  of  the  three  fluid 
media,  respectively.  Both  the  ocean  surface  and  ocean  bottom  are  modeled  as  plane,  parallel 
boundaries,  and  die  ocean  is  D  meters  deep.  The  source  distribution  is  a  motionless,  time- 
harmonic,  point  source  [see  (2.1-8)  and  (2.1-11)]  located  in  medium  II  at  horizontal  range  r  =  0 
meters  and  depth  y  =  y0  meters  where  the  cylindrical  coordinates  (r,  (p,y)  are  illustrated  in  Fig.  2.1- 
3. 

If  the  ocean  surface  is  modeled  as  an  ideal  pressure-release  boundary  and  the  ocean  bottom 
is  modeled  as  an  ideal  rigid  boundary,  then  it  can  be  shown  that  the  time-invariant,  space-variant, 
complex  frequency  response  of  the  ocean  at  frequency  /  hertz  is  given  by  [6] 


(2.1-33) 


where 


Np  =  INTI 


4  Df 


l 

V2L  v2 


—  1 


+  1 


(2.1-34) 


is  the  total  number  of  propagating  normal  modes  excited  by  the  source  frequency  /  hertz, 
is  the  zeroth-order  Hankel  function  of  the  second  kind. 


f<fn . 


(2.1-35) 


is  the  propagation-vector  component  in  the  horizontal  radial  direction,  in  radians  per  meter,  of  the 
nth  normal  mode, 


k2  =  2nf  I  c2  =  2  njX2  (2.1-36) 

is  the  wavenumber  in  medium  II  in  radians  per  meter, 
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I  p,c, 


Figure  2.1-2.  Simple  waveguide  model  of  the  ocean  where  pxcv  p2c2,  and  p3c3  are  the 
characteristic  impedances  of  fluid  media  I  (air),  II  (seawater),  and  III  (ocean  bottom), 
respectively.  The  source  distribution  is  a  motionless,  time-harmonic,  point  source  located  in 
medium  II  at  horizontal  range  r  =  0  meters  and  depth  y  =  y0  meters. 


10 


(2n  +  l)c2 
Jn  4D 


(2.1-37) 


is  the  cutoff  frequency  in  hertz  of  the  nth  normal  mode,  and 


k 


(2n  +  \)n 
2D 


(2.1-38) 


is  the  propagation- vector  component  in  the  depth  direction,  in  radians  per  meter,  of  the  nth  normal 
mode.  Note  that  in  the  far-field,  that  is,  when  K2nr>>  1, 


exp 


r 


-j 


i  71 

kr  r - 

2-"  4 


kr  r»  1 . 

r2.n 


(2.1-39) 


By  inspecting  (2.1-33),  it  can  be  seen  that  the  complex  frequency  response  depends  on  the  actual 
values  of  the  depth  coordinates  of  the  field  point  (y)  and  the  point  source  (y0)  and  not  on  their 
difference.  This  is  an  indication  that  the  medium  is  space-variant  in  the  depth  direction. 

A  second  example  of  a  time-invariant,  space-variant,  complex  frequency  response  can  be 
obtained  by  reconsidering  the  simple  waveguide  model  of  the  ocean  illustrated  in  Fig.  2.1-2.  If  the 
ocean  surface  is  once  again  modeled  as  an  ideal  pressure-release  boundary  as  before,  and  if  the 
ocean  bottom  is  now  modeled  as  the  boundary  between  two  different  fluid  media  instead  of  being 
modeled  as  an  ideal  rigid  boundary,  then  it  can  be  shown  that  the  time-invariant,  space-variant, 
complex  frequency  response  of  the  ocean  at  frequency  /  hertz  for  a fast  fluid-bottom  is  given  by 
[7] 


Hu( /.r|r„)  =  Hu{f,r,*,y |y„)  =  -jj £  E;' 0<y<D, 


(2.1-40) 


where 


f ,  r 


N,  =  INTI 


V~L 


4Dfcos/3c 


-1 


+  1 


(2.1-41) 


is  the  total  number  of  trapped  (propagating)  normal  modes  excited  by  the  source  frequency  /  hertz, 

ffc  =  sin-1  ( c2  /c3 ) ,  c3  >  c2 ,  (2.1  -42) 

is  the  critical  angle  of  incidence, 

2 ky  D  -  sin(2 £>)  -  2  ( p2 /p3 )3  tan(^  £>)sin2(^ D) 


E  = 


4kv 


(2.1-43) 
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is  the  “energy”  in  meters  of  the  nth  eigenfunction,  Hq2\»)  is  the  zeroth-order  Hankel  function  of 
the  second  kind, 


kr2,a  =  k2  sin&„  (2.1-44) 

is  the  propagation-vector  component  in  the  horizontal  radial  direction,  in  radians  per  meter,  of  the 
nth  normal  mode  (see  Fig.  2.1-4), 


ky2  *  -  k2  cos  fl2n  (2.1-45) 

is  the  propagation- vector  component  in  the  depth  direction,  in  radians  per  meter,  of  the  nth  normal 
mode  (see  Fig.  2.1-4), 


k2  =  2nf  / c2  =  2n/X2 


(2.1-46) 


is  the  wavenumber  in  medium  II  in  radians  per  meter,  and  the  angles  of  propagation  fl2n , 
n  =  0,1,...,  Nt  - 1,  (see  Fig.  2.1-4)  are  the  roots  of  the  following  transcendental  equation: 


tan 


2  nf 

V  C2 


D  cos  fi2n 


p2c3  cos(32n 


p2c2A/(sin/32n  /sin/3c)2  - : 


=  0,  c3>c2,  f3c  <  P2  „  <n/2. 


(2.1-47) 


The  far-field  approximation  of  Hq2)(*)  is  given  by  (2.1-39).  By  inspecting  (2.1-40),  it  can  be  seen 
that  the  complex  frequency  response  depends  on  the  actual  values  of  the  depth  coordinates  of  the 
field  point  (y)  and  the  point  source  (y0)  and  not  on  their  difference.  This  is  an  indication  that  the 
medium  is  space-variant  in  the  depth  direction. 

A  third  and  final  example  of  a  time-invariant,  space-variant,  complex  frequency  response  can 
be  obtained  by  reconsidering  the  simple  waveguide  model  of  the  ocean  illustrated  in  Fig.  2.1-2.  If 
both  the  ocean  surface  and  the  ocean  bottom  are  modeled  as  boundaries  between  two  different  fluid 
media,  then  it  can  be  shown  that  the  time-invariant,  space-variant,  complex  frequency  response  of 
the  ocean  at  frequency  /  hertz  is  given  by  [8] 


HM(f,r\r0)=Hu(f,r,<t>,y\y0) 

.  1 

=  +J~r~ 

(  r»y0)  jfl2i  expf-^j,  (y  +  y0)]  +  exp[+^2  (y -  %)]} 
o  *r2 

xj0(krr)krdkr,  0<y<y0, 

(2.1-48) 

when  the  field  point  is  at  a  depth  between  the  ocean  surface  and  the  point  source  ( 0  <  y  <  y0),  and 
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Figure  2.1-4.  The  propagation  vector  k2„  of  the  nth  trapped  (propagating)  normal  mode  in 
medium  II.  Note  that  |k2  „|  =  k2  =  2; vf/c2 .  Also  shown  are  the  angle  of  propagation  fi2n  and  the 
scalar  propagation-vector  components  kr  and  kY  ^  in  the  horizontal  range  r  and  depth  Y 
directions,  respectively. 


HM{f,r\r0)=  HM(f,r,<l),y\y0) 

.  1 

-  +J~r~ 
4nJ 

\r'y) {p2i exp[  JK (y + %)] + exp[-j\ (y - y0)]} 

0  * y2 

xj0(krr)krdkr,  y0<y<D, 

(2.1-49) 

when  the  field  point  is  at  a  depth  between  the  point  source  and  the  ocean  bottom  (y0<y<D), 
where 


F(kr,y)  = 


1  +  fl23  expj-y'2^  (D  -  y) 
l-R21R23exp(-j2kYD)  ’ 


(2.1-50) 


and 


F(K.y«)  =  nK,y)\r.n, 
r  —  Ei  ~  El 

P\ky2  +  Pi^ 


(2.1-51) 

(2.1-52) 


K  _  PiKi  Ei-h 
Pb  ^r,  4  P2  ky. 


(2.1-53) 


are  the  plane-wave  velocity  potential  reflection  coefficients  at  the  ocean  surface  and  ocean  bottom, 
respectively. 


V*?-*?. 

k2r  <k]. 

i  =  1,2,3, 

*?>*?. 

i  =  1,2,3, 

(2.1-54) 


are  the  propagation-vector  components  in  the  depth  direction,  in  radians  per  meter,  in  media  I,  D, 
and  HI,  respectively, 

k i  =  2nf/ci ,  i  =  1,2,3,  (2.1-55) 


are  the  wavenumbers  in  radians  per  meter  in  media  1, 13,  and  HI,  respectively,  and  J0(«)  is  the 
zeroth-order  Bessel  function  of  the  first  kind.  Note  that  in  the  far-field,  that  is,  when  krr»  1 , 


J0(Kr)~ 


cos 


r 


krr»  1. 


(2.1-56) 
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Note  that  the  function  F(kr,y)  in  the  integrand  of  (2.1-49)  is  not  evaluated  at  y  =  y0.  By  inspecting 
(2.1-48)  and  (2.1-49),  it  can  be  seen  that  the  complex  frequency  response  depends  on  the  actual 
values  of  the  depth  coordinates  of  the  field  point  (y)  and  the  point  source  (y0)  and  not  on  their 
difference.  This  is  an  indication  that  the  medium  is  space-variant  in  the  depth  direction. 


Example  2.1-3  Time-Variant,  Space-Invariant,  Ocean-Medium  Filters 

If  the  ocean-medium  filter  is  time-variant  and  space-invariant,  then  [1] 

fytf(^’rko>ro)  =  — *0  |^o)> 


(2.1-57) 


that  is,  the  impulse  response  depends  on  the  vector  spatial  difference  R  =  r  -  r0.  The  vector  spatial 
difference  R  corresponds  to  the  distance  and  direction  between  source  and  field  points.  The 
magnitude  of  the  vector  spatial  difference  |R|  =  |r  -  r0|  is  a  measure  of  how  far  away  the  impulse 
was  applied,  or  the  distance  traveled  between  source  and  field  points.  Substituting  (2.1-57)  into 
(2.1-2)  yields 

y M  (*> r)  =  XM  (to  ’  r0  (t’ r  —  ro  I  )dto  dr, . 

In  addition,  substituting  (2.1-57)  into  (2.1-4)  yields 

Hu(t,r\f,r0)  =  j"mhM(t,r-  r0 1  t0)exp[-j2nf(t  -  t0)]dt0 
=  Hu(t,r-r0\f). 


(2.1-58) 


(2.1-59) 


If  the  source  distribution  is  motionless  and  time-harmonic  [see  (2.1-8)]  and  the  ocean- 
medium  filter  is  time-variant  and  space -invariant,  then  substituting  (2.1-59)  into  (2.1-10)  yields  the 
following  expression  for  the  output  acoustic  signal  (velocity  potential)  from  the  ocean-medium 
filter: 


yM(t, r)  =  x/, « (O # « r  -  £ | /) dZ,txp(+j2nft) 

=  [^,  M  (r)  *  Hm  (t,  r  |  /)  j  cxp(+j2nft) , 


(2.1-60) 


where  the  time-variant,  space-invariant,  complex  frequency  response  of  the  ocean  at  frequency  / 
hertz  is  given  by  [see  (2.1-59)] 


HM(t,  r  -  r0 1/)  =  \~JiM(t,r-  r0 1  t0)exp[-j2nf(t  -  t0)]dt0. 


(2.1-61) 
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Furthermore,  if  the  source  distribution  is  a  motionless,  time-harmonic,  point  source,  then  (2.1-8) 
and  (2.1-11)  are  applicable,  and  upon  substituting  either  (2.1-11)  into  (2.1-60),  or  (2.1-59)  into 
(2.1-12),  we  obtain 


3W(*,r)  =  S0HM(t,  r-r0  \f)exp(+j2nft), 


(2.1-62) 


where  HM(t, r-r0 1/)  is  given  by  (2.1-61). 


2.2  Pulse  Propagation  In  The  Ocean 

Before  we  can  discuss  the  pulse-propagation  problem  in  this  section,  we  need  to  introduce 
two  coupling  equations.  The  first  coupling  equation  is  the  transmit-aperture  coupling  equation  and 
is  given  by  (see  [9]  and  Fig.  2.2-1) 


xM(t, r)  =  f  X(f,r)AT(f,r)txp(+j2nft)df, 
•1—00 


(2.2-1) 


where  xM(t,r)  is  the  input  acoustic  signal  (source  distribution)  to  the  fluid  medium  with  units  of 
inverse  seconds  (volume  flow  rate  per  unit  volume  of  fluid)  at  time  t  and  position  r  =  (x,y,z)  of  the 
transmit  aperture, 


X(f,r)  =  Ft  {*(?,  r)}  =  J__  x(t,  r)  exp(-j2nft)  dt  (2.2-2) 


is  the  complex  frequency  spectrum  of  the  input  electrical  signal  at  position  r  =  (x,y,z)  of  the 
transmit  aperture  with  units  of  volts  per  hertz,  x(t,r)  is  the  input  electrical  signal  to  the  transmit 
aperture  at  time  t  and  position  r  =  (x,y,z)  with  units  of  volts,  and  AT(f,  r)  is  the  complex 
frequency  response  (amplitude  and  phase)  at  position  r  =  (x,y,z)  of  the  transmit  aperture,  also 
known  as  the  complex  transmit  aperture  function,  with  units  of  inverse-seconds  per  volt,  where  f 
represents  input  or  transmitted  frequencies  in  hertz.  Equation  (2.2-1)  represents  the  coupling  of  the 
input  electrical  signal  to  the  fluid  medium,  that  is,  the  production  of  the  input  acoustic  signal  (source 
distribution)  due  to  the  input  electrical  signal  via  the  transmit  aperture  whose  performance  is 
described,  in  part,  by  the  complex  transmit  aperture  function.  Another  useful  result  can  be  obtained 
from  (2.2-1)  as  follows.  Since 


xM(t,r)  =  Ff-l{X(f,r)AT(f,r)}, 


(2.2-3) 


then 


XM(f, r)  =  /){*M(?,r)}  =  j  ~xM(t,r)exp(-j2nft)dt 
=  X(f,r)AT(f,r), 


(2.2-4) 


or. 
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x(t,r) 


X(f,r) 


Ar(f,r) 

/jw(t,r|r0,r0) 

AR(ri,  r) 

W r) 

YM(V,rf 

Y«l,  r) " 

TRANSMIT 

APERTURE 

MEDIUM 

FILTER 

RECEIVE 

APERTURE 

Figure  2.2-1.  Block  diagram  model  of  small-amplitude  pulse  propagation. 
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XM(f,r)  =  X(f,r)AT(f,r), 


(2.2-5) 


where  XM(f,  r)  is  the  complex  frequency  spectrum  of  the  input  acoustic  signal  (source 
distribution)  at  position  r  =  (x,  y,z)  of  the  transmit  aperture  with  units  of  inverse-seconds  per  hertz 
(see  [9]  and  Fig.  2.2-1). 

The  second  coupling  equation  is  the  receive-aperture  coupling  equation  and  is  given  by  (see 
[9]  and  Fig.  2.2-1) 

y(r,r)  =  P  YM(ri,r)AR{r),r)exp{+j2nr}t)dr],  (2.2-6) 

J  —oo 

where  y(t,  r)  is  the  output  electrical  signal  from  the  receive  aperture  at  time  t  and  position 
r  =  ( x,y,z )  with  units  of  volts, 

=  i;{yw(r,r)}  =  J~  yM(t,r)exp(-j2nr]t)dt  (2.2-7) 

is  the  complex  frequency  spectrum  of  the  acoustic  signal  (velocity  potential)  incident  upon  the 
receive  aperture  at  position  r  =  (x,y,z)  with  units  of  (m2/sec)/Hz,  yM(t, r)  is  the  acoustic  signal 
(velocity  potential)  incident  upon  the  receive  aperture  at  time  t  and  position  r  =  (x,y,  z)  with  units  of 
squared-meters  per  second,  and  AR(ri,r )  is  the  complex  frequency  response  (amplitude  and  phase) 
at  position  r  =  (x,y,z)  of  the  receive  aperture,  also  known  as  the  complex  receive  aperture  Junction, 

with  units  of  V/(m2/sec),  where  rj  represents  output  ox  received  frequencies  in  hertz.  Note  that  in 

general,  Equation  (2.2-6)  represents  the  coupling  of  the  fluid  medium  to  the  output  electrical 

signal,  that  is,  the  production  of  the  output  electrical  signal  due  to  the  input  acoustic  signal  via  the 
receive  aperture  whose  performance  is  described,  in  part,  by  the  complex  receive  aperture  function. 
Another  useful  result  can  be  obtained  from  (2.2-6)  as  follows.  Since 

y(t,r)  =  F,"‘{l'„(>I,r)AJ(l),r)},  (2.2-8) 

then 

Y(ri,r)  =  ^{y(f,r)}  =  jj(t,r)exp(-j2ni]t)dt 
=  YM(r],r)AR(q,r), 

or, 

Y(r],  r)  =  YM(T!,r)  AR(r],r),  (2.2  - 10) 

where  Y(r), r)  is  the  complex  frequency  spectrum  of  the  output  electrical  signal  at  position 
r  =  (x,  y,  z )  of  the  receive  aperture  with  units  of  volts  per  hertz  (see  [9]  and  Fig.  2.2-1). 
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Now  that  we  have  introduced  the  two  coupling  equations,  we  will  next  show  how  to  use  them 
to  solve  the  pulse-propagation  problem.  Recall  from  Section  2.1  that  xM(r,r)  and  yw(?,r)  are 
related  by  the  following  linear,  three-dimensional,  inhomogeneous  wave  equation  [see  (2.1-1)]: 

V2yM(?,r)-^^|jyw(t,r)  =  xM(t,  r).  (2.2-11) 

Also  recall  from  Section  2.1  that  the  solution  of  (2.2-11)  is  given  by  [see  (2.1-2)] 

=  /  M  j"  „ ^ ’ r0 (^’ I  ’ *0 )dt0 dr0 ,  (2.2-12) 

where  /zA/(t,r  |r0,r0)  is  the  time-variant,  space-variant,  impulse  response  (Green’ s  function)  of  the 
medium  filter  (see  Fig.  2.2-1),  that  is,  it  is  the  response  of  the  filter  at  time  t  and  position 
r  =  (x,y,z)  due  to  the  application  of  a  unit-amplitude  impulse  at  time  t0  and  position 
r0  =  (x0,y0,z0).  Substituting  the  first  coupling  equation  given  by  (2.2-1)  into  (2.2-12)  and  making 
use  of  (2.1-4)  yields 


yM(t,r)  =  J_~ f_~  X(f,r0)AT(f,r0)HM(t,r\f,r0)  dr0  zxp(+j2nft)df. 


(2.2-13) 


where 


Hm  (r,  r  |  /, r0 )  =  J“  hM  (t, r  1 10 , r0 )  &xp[-j2nf(t  - 10 )] dt0  (2.2- 14) 

is  the  time-variant,  space-variant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz, 
[see  (2.1-4)  and  Fig.  2.2-1],  Note  that  since  the  impulse  response  of  the  ocean  medium 

hM (t, r|r0,r0)  is  a  real  function, 


I  f^0)-HM(t,r\f,r0). 


(2.2-15) 


Equation  (2.2-13)  is  a  general  solution  for  the  acoustic  signal  (velocity  potential)  that  is 
incident  upon  the  receive  aperture.  It  can  be  used  in  conjunction  with  (2.2-6)  and  (2.2-7)  to  solve 
for  the  output  electrical  signal  from  the  receive  aperture.  Before  we  discuss  how  to  simplify  (2.2- 
13)  in  order  to  solve  real-world,  pulse-propagation  problems,  let  us  make  the  following  interesting 
observation.  If  we  refer  back  to  the  time-harmonic  solution  given  by  (2.1-10),  and  if  we  set 
xf  M( C)  =  X(fX)AT(f,Q,  then  (2.2-13)  can  be  interpreted  as  being  equal  to  the  integration  of  the 

time-harmonic  solution  given  by  (2.1-10)  over  all  the  frequency  components  contained  in  the 
complex  frequency  spectrum  of  the  input  electrical  signal  applied  to  the  transmit  aperture. 

We  begin  the  simplification  of  (2.2-13)  by  assuming  that  an  identical  input  electrical  signal  is 
applied  at  all  spatial  locations  r0  of  the  transmit  aperture,  that  is, 

x(r,r0)  =  x(0  (2.2-16) 


so  that  [see  (2.2-2)] 
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X(f,r0)-X(f). 


(2.2-17) 


Equations  (2.2-16)  and  (2.2-17)  are,  in  fact,  valid  in  the  vast  majority  of  real-world,  practical 
problems.  For  example,  in  most  practical  situations,  an  identical  input  electrical  signal  is  applied  to 
all  elements  in  a  transmit  array  before  the  application  of  complex  weights  for  beamforming 
purposes.  Substituting  (2.2-17)  into  (2.2-13)  yields 


yw(*,r)  =  j_~X(f)j_~AT(f,r0)HM(t,r\f,r0)dr0exip(+j27tft)df.  (2.2-18) 

Equation  (2.2-18)  can  be  simplified  by  specifying  the  transmit  aperture.  In  this  paper,  both  the 
transmit  and  receive  apertures  will  be  modeled  as  volume,  conformal  arrays  composed  of  unevenly- 
spaced,  complex-weighted,  point  elements.  This  type  of  model  for  both  of  the  apertures  will  give  us 
maximum  flexibility.  Therefore,  the  complex  transmit  aperture  function  is  given  by 

Ar(/,r0)  =  %Ti  ( f)STi  (/)5(  r0  -  ^ ),  (2.2-19) 

i=i 

where  NT  is  the  total  number  of  elements  in  the  transmit  array, 

cr,  (/)  =  aTXf)exp[+jdTif)]  (2.2-20) 

is  the  dimensionless,  frequency-dependent  complex  weight  associated  with  element  i  in  the  transmit 
array,  where  aT(f)  and  0TXf)  are  real,  frequency-dependent,  amplitude  and  phase  weights, 

respectively,  STj  (f)  is  the  frequency-dependent  transmitter  sensitivity  of  element  i  in  the  transmit 

array  with  units  of  (m3/ sec) /V  ,rT.  =  (xTj  ,yT.,zTj)  is  the  position  vector  to  element  i  in  the  transmit 
array,  and  the  impulse  function  <5(r0  -iy),  with  units  of  inverse  cubic  meters,  represents  a  point 
element  at  rr  =  (xT  ,yT  ,zT).  The  real  amplitude  weights  are  dimensionless  and  are  even  functions 
of  frequency,  that  is, 


aTX-f)  =  aTif), 


(2.2-21) 


and  the  real  phase  weights,  with  units  of  radians,  are  odd  functions  of  frequency,  that  is, 


eTx-f)=-eTxf). 

Therefore, 


c'TX-f)  =  cTXf). 


(2.2-22) 


(2.2-23) 


In  addition,  since  the  point  elements  are  electroacoustic  transducers  with  real  impulse-response 
functions, 


s;x-f)=sTxf). 


(2.2-24) 
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Because  of  (2.2-23)  and  (2.2-24), 


A*T(-f,r0)  =  AT(f,r0). 


(2.2-25) 


Substituting  (2.2-19)  into  (2.2-18)  yields 


y„  it,  r)  =  J~  X(  f )  £  cTj  (  / )  ST ( f )  Hm  ( t,  r  |  /,  rTt )  exp(+j  In  ft)  df ,  (2.2-26) 

1  =  1 

or 

yM(t,r)  =  E  12  ^(  /  )cr,  ( / )  ^7-  ( / )  r  I  /’  rr, )  exp(+ j2^//)  df .  (2.2-27) 

1=1 

Equation  (2.2-27)  indicates  that  the  acoustic  signal  (velocity  potential)  yM(t, r)  incident  upon  the 
receive  aperture  at  time  t  and  position  r  =  (x,y,z)  is  equal  to  the  sum  of  the  contributions  from  each 
element  in  the  transmit  array. 

Equation  (2.2-27)  can  be  simplified  by  specifying  the  class  of  transmitted  electrical  signals  to 
be  considered.  Typical  transmitted  electrical  signals  used  in  active  radar  and  sonar  systems  can  be 
modeled  as  amplitude-and-angle-modulated  carriers  as  follows  [10]: 

x(t)  =  Aa{t)cos[2n  fct  +  0(O]rect^- — (2.2-28) 

where  A  is  the  amplitude  scaling  factor  with  units  of  volts,  a(t)  is  a  dimensionless  amplitude- 
modulating  signal  where  |a(t)|  <  1,  fc  is  the  carrier  frequency  in  hertz,  6(t)  is  an  angle-modulating 
signal  (also  known  as  the  phase  deviation)  with  units  of  radians,  cos(2n ft)  is  the  carrier 
waveform,  and 


rect 


ft-o.5r 

V  T  J 


0  <t<T, 
otherwise, 


(2.2-29) 


is  the  rectangle  function  where  T  is  the  pulse  length  in  seconds.  The  corresponding  complex 
envelope  of  x(t),  denoted  by  x(t),  is  given  by  [10] 


x(t)  -  A  a(t)  exp[+ j0(t)]  rect 


ft_  -0.57^ 

t  y 


(2.2-30) 


provided  that  f  >  W,  where  W  is  the  bandwidth  of  x(t)  in  hertz.  The  amplitude-and-angle- 
modulated  carrier  x(t)  is  a  real,  bandpass  signal  whose  magnitude  spectrum  |X(/)|  is  centered 
about  /  =  ±  f  hertz,  where 


X(f)  =  7^{^(0}  =  j_^x(t)exp(-j27tft)dt .  (2.2-31) 
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However,  its  complex  envelope  x(t)  is  a  complex ,  baseband  signal  whose  magnitude  spectrum 
|x(/)|  is  centered  about  /  =  0  hertz,  where 


X(f)  =  *?{*(*)}  =  [j(t)exp(-j27tft)dt.  (2.2-32) 


The  relationship  between  X(  f)  and  X(f)  is  given  by  [1 1] 


X(/)  =  0.5[*(/-/c) +  *•(-[/ +  /J)]. 


(2.2-33) 


Substituting  (2.2-33)  into  (2.2-27)  and  making  use  of  (2.2-15),  (2.2-23),  and  (2.2-24)  yields 


yM(t,r)  =  Re{  yM  (t,  r)  exp(+j2nfct)} ,  (2.2-34) 


where 


|r.  X{f)cTj(f  +  fc)ST(f  +  fc)HM[t,r\f  +  fc,rT^exp(+j27tft)df 

(2.2-35) 

is  the  complex  envelope  of  yM(t,r).  Equations  (2.2-34)  and  (2.2-35)  describe  the  acoustic  signal 
(velocity  potential)  yM(t,r)  that  is  incident  upon  the  receive  aperture  at  time  t  and  position 
r  =  (x,y,z)  in  terms  of  its  complex  envelope  yM(t,r).  What  we  need  to  do  next  is  to  turn  our 
attention  to  the  receive  aperture  itself  so  that  we  can  predict  the  output  electrical  signal  from  the 
receive  aperture. 

Recall  that  both  the  transmit  and  receive  apertures  are  being  modeled  as  volume,  conformal 
arrays  composed  of  unevenly-spaced,  complex-weighted,  point  elements.  Therefore,  the  complex 
receive  aperture  function  is  given  by 


ar  (v, r)  =  X  < CR„  07)  ■ SR„  07)  < £( : r  “  rRn )  ■ .  (2.2-36) 

n=l 

where  NR  is  the  total  number  of  elements  in  the  receive  array, 

cr„  (V)  =  aRn  (r?)exp  [+j0Rt  (??)]  (2.2-37) 

is  the  dimensionless,  frequency-dependent  complex  weight  associated  with  element  n  in  the  receive 
array,  where  aRfrj)  and  dR  (V)  are  real,  frequency-dependent,  amplitude  and  phase  weights, 

respectively,  SR  (77)  is  the  frequency-dependent  receiver  sensitivity  of  element  n  in  the  receive 

array  with  units  of  (V-m3)/(m2/sec),  or  Vy/((m2/sec)/m3),  rRn  =  {xR^,yR^,zR)  is  the  position 
vector  to  element  n  in  the  receive  array,  and  the  impulse  function  S(r-rR  ),  with  units  of  inverse 
cubic  meters,  represents  a  point  element  at  rR„  =  (xR^,yR^,zRn )■  The  real  amplitude  weights  are 
dimensionless  and  are  even  functions  of  frequency,  that  is, 
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or 

yit,rK' )  =  ^_YM(r],rR  )cR  ^)SR  (t])exp(+j2nrit)drj,  n  =  1,2 ,...,NR, 

(2.2-44) 

where  y(t,rR  )  is  the  output  electrical  signal  from  element  n  in  the  receive  array  at  time  t  and 
position  )  with  units  of  volts, 

(2.2-45) 

is  the  complex  frequency  spectrum  of  the  acoustic  signal  (velocity  potential)  incident  upon  element 
n  in  the  receive  array  at  7]  hertz  and  position  rR„  =  (xRa,yR,zR)  with  units  of  (m2/sec)/  Hz  [see 
(2.2-7)], 
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yM(t,rRn )  =  Re{  yM  (t,  )  exp(+;2^/c  *)} 


(2.2-46) 


is  the  acoustic  signal  (velocity  potential)  incident  upon  element  n  in  the  receive  array  at  time  t  and 
position  rR __  =  (xR„  ,yRn,zRn )  with  units  of  squared-meters  per  second  [see  (2.2-34)],  and 


y«  «.■*.)= X  £  ■ *( . / )  < h  ( / + f, )  ■ s,  ( /  ■ +  /, ) « *  1 
1=1 

f  +  fc^T, 

exp{+jl7tft)df 

(2.2-47) 

is  the  complex  envelope  of  yM(t, r)  [see  (2.2-35)].  Equations  (2.2-44)  through  (2.2-47)  are  the 
pulse-propagation  coupling  equations.  They  couple  the  output  electrical  signal  at  a  point  element  in 
a  receive  array  to  the  transmitted  electrical  signal  at  the  input  to  a  transmit  array  via  the  complex 
frequency  response  of  a  fluid  medium  (e.g.,  air  or  water).  The  pulse-propagation  coupling 
equations  are  based  on  linear,  time-variant,  space-variant,  filter  theory,  the  principles  of  complex 
aperture  theory  and  array  theory,  and  solving  the  linear  wave  equation.  They  provide  a  consistent, 
logical,  and  straightforward  mathematical  framework  for  the  solution  of  small-amplitude,  acoustic, 
pulse-propagation  problems.  Examples  on  the  use  of  the  pulse-propagation  coupling  equations  are 
given  in  Section  3. 


3  Bistatic  Scattering 
3.1  No  Motion 

In  this  section  we  will  analyze  the  bistatic  scattering  problem  shown  in  Fig.  3.1-1.  The 
transmitter,  the  discrete  point  scatterer,  and  the  receiver  are  not  in  motion.  For  example  purposes, 
they  are  placed  in  deep  water  so  that  the  ocean  medium  can  be  treated  as  being  unbounded,  that  is, 
there  are  no  boundaries  and,  hence,  no  reflections.  The  speed  of  sound  and  ambient  density  are 
constants.  Although  the  propagation  of  sound  between  the  source  and  discrete  point  scatterer,  and 
between  the  discrete  point  scatterer  and  receiver  can  be  treated  as  transmission  through  linear,  time- 
invariant,  space-invariant  filters;  the  overall  solution  for  this  bistatic  scattering  problem  corresponds 
to  transmission  through  a  linear,  time-invariant,  space-variant  filter.  The  presence  of  a  discrete  point 
scatterer  in  an  unbounded,  homogeneous  fluid  medium  (i.e.,  a  fluid  medium  with  constant  speed  of 
sound  and  ambient  density)  causes  the  medium  to  be  space- variant. 

Let  the  source  distribution  xM(t,  r)  be  a  motionless,  time-harmonic,  point  source  with  units  of 
inverse  seconds,  that  is,  let  [see  (2.1-8)  and  (2.1-11)] 

xM(t,  r)  =  S0S(r  -  r0)exp(+j2nft),  (3.1-1) 

where  S0  is  the  source  strength  in  cubic  meters  per  second  and  the  impulse  function  <5(r  -  r0),  with 
units  of  inverse  cubic  meters,  represents  a  point  source  at  r0  =  (x0,y0,z0).  The  propagation  of 
sound  between  the  source  and  discrete  point  scatterer  can  be  modeled  as  transmission  through  a 
linear,  time-invariant,  space-invariant,  filter.  Therefore,  the  acoustic  field  (velocity  potential) 
incident  upon  the  discrete  point  scatterer  at  r,  =  (x^y^z,)  is  given  by  [see  (2.1-17),  (2.1-20),  and 
(2.1-22)] 
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y*(Mi)  =  y/tj#(ii)«p(+72jr/0. 


(3.1-2) 


where 


expf-yfclr.-rJ) 

yf.M  Or.)  =  ^0  -  r0)  =  -  ^  -  •  (3-1_3) 

47r|1i-ro| 

As  can  be  seen  from  Fig.  3.1-1,  the  position  vector  from  the  point  source  to  the  discrete  point 
scatterer  is  given  by 


ro,i  =  1i-ro- 


Therefore,  (3.1-3)  can  be  rewritten  as 


)V,M(ri)  —  ~  S0 


exp| 

M' 

r°.'l) 

4  n 

ro. 

l] 

(3.1-4) 


(3.1-5) 


In  order  to  compute  the  acoustic  signal  incident  upon  the  receiver,  we  treat  the  discrete  point 
scatterer  as  another  motionless,  time-harmonic,  point  source  with  units  of  inverse  seconds,  that  is, 
let  [see  (3.1-1)  and  Fig.  3.1-2] 

x'M(t,  r)  =  SqS(t  -  r,  )exp(+j'27r/0,  (3.1-6) 

where  S'0  is  the  source  strength  in  cubic  meters  per  second  and  the  impulse  function  <5(r  -  Tj) ,  with 
units  of  inverse  cubic  meters,  represents  a  point  source  at  Tj  =  (Xp^pZ,).  The  source  strength  is 
given  by 


where  yf  M(rl)  is  the  spatial-dependent  part  of  the  time-harmonic  velocity  potential  at  r,  with  units 

of  squared-meters  per  second,  given  by  either  (3.1-3)  or  (3.1-5),  and  g,( /^o.i’^u)  is  the 
scattering  amplitude  function  of  the  discrete  point  scatterer  with  units  of  meters.  The  unit  of  meters 
for  the  scattering  amplitude  function  represents  an  effective  scattering  length  that  may  be  larger  or 
smaller  than  the  actual  length  of  the  scatterer.  The  scattering  amplitude  function  is  a  complex 
function  (magnitude  and  phase)  and  is,  in  general,  a  function  of  frequency,  the  direction  of  wave 
propagation  from  the  source  to  the  scatterer  (n01),  and  the  direction  of  wave  propagation  from  the 
scatterer  to  the  receiver  (n,  2 )  (see  either  Fig.  3.1-1  or  Fig.  3.1-2)  [12]. 

The  propagation  of  sound  between  the  discrete  point  scatterer  and  receiver  can  also  be 
modeled  as  transmission  through  a  linear,  time-invariant,  space-invariant  filter.  Therefore,  the 
acoustic  field  (velocity  potential)  incident  upon  the  receiver  at  r2  =  (x2,y2,z2)  due  to  a  point  source 
at  q  =  (x,,y,,z,)  is  given  by  [see  (2.1-17),  (2.1-20),  and  (2.1-22)] 

3;M(^r2)  =  3;/,A/(r2)exP(+727z:/?),  (3.1-8) 


where 


S'<S(r-ri) 


Figure  3.1-2.  Bistatic  scattering  geometry.  Both  the  transmitter  at  point  0,  P0(r0),  and  the 
discrete  point  scatterer  at  point  1,  are  time-harmonic  point  sources.  None  of  the  platforms 

are  in  motion. 


(3.1-9) 


v  fri  -yrr  tfr  ,1  „«p(-.i*l1i-1il) 

yf,M\r2)~\^M\J'<r2  ri)~  4^:|r  —  r| 


As  can  be  seen  from  either  Fig.  3.1-1  or  Fig.  3.1-2,  the  position  vector  from  the  discrete  point 
scatterer  to  the  receiver  is  given  by 


ri,2  =  r2_1l- 


Therefore,  (3.1-9)  can  be  rewritten  as 


expf- jk\r.  2 1 ) 

yf.M(h)  =  -WM,2)  =  -S'  , -  t . * 

47rru 


(3.1-10) 


(3.1-11) 


Let  us  now  begin  the  process  of  obtaining  final  expressions  for  both  the  time-harmonic 
velocity  potential  and  acoustic  pressure  incident  upon  the  receiver.  Substituting  (3.1-7)  and  (3.1-5) 
into  (3.1-11)  yields 


yf,M(r2)  =  S0HM(f^0,l)8l(fAlA2)HM(f^U)’ 


or,  equivalently. 


where 


y/,M^r2^  —  S0Hm(  f,r2 1 r0), 

( /.  r,  |  r0 )  =  Hm  ( /,  r0>1 )  g,  ( /,  n0>1 ,  nu )  Hu  ( /,  r12 ) 


exp 

~jk( 

ro,i 

+ 

ru|). 

167T2 

ro,i 

ru 

>  1 

(  .  \  exP[-7^(|ri  ~  ro| +  |r2  “  ri|)] 

=  gi  />«0,P«1,2  - 7731 -  r  „  „ - ■ 

16^  |ri-r0  r2_ri 


(3.1-12) 

(3.1-13) 

(3.1-14) 


Note  that  if  the  bis  tan  c  scattering  problem  shown  in  Fig.  3.1-1  corresponded  to  transmission 
through  a  space-invanant  filter,  then  the  complex  frequency  response  given  by  (3.1-14)  would  be  a 
function  of  the  vector  spatial  difference  r2  -  r0,  which  it  is  not.  Also  note  that  in  order  to  model  the 

effects  of  frequency-dependent  attenuation,  simply  replace  the  real  wavenumber  k  given  by  (2.1- 
23)  with  the  following  complex  wavenumber  K : 


K  =  k- ja(f). 


(3.1-15) 


where  Ct(f)  is  the  real,  frequency-dependent,  attenuation  coefficient  in  nepers  per  meter.  In 
addition  to  being  real  quantities,  both  k  and  a(f )  are  positive. 
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With  the  use  of  (2.1-23),  (3.1-8),  (3.1-13),  and  (3.1-14);  and  by  replacing  the  real 
wavenumber  k  in  (3.1-14)  with  die  complex  wavenumber  K  given  by  (3.1-15),  we  can  summarize 
our  results  as  follows:  for  the  bistatic  scattering  problem  shown  in  Fig.  3.1-1,  the  time-harmonic 
velocity  potential  in  squared-meters  per  second  incident  upon  the  receiver  at  r2  =  (x2,y2,z2)>  due  to 
a  time-harmonic  point  source  at  r0  =  (x0,y0,z0)  and  a  discrete  point  scatterer  at  r,  =  (xi,y1,zi),  is 
given  by 


yM(t,r2)  =  S0 Hm ( /, r2 1  r0 ) exp(+j2nft),  t>  r , 


(3.1-16) 


where  S0  is  the  source  strength  in  cubic  meters  per  second, 


,  ,  ,  ,  .  „  ,  exp[-a(/)(r0l  +  ru)l 

Hu(f,t2  =  ,  - iexp(-;2^/r) 


(3.1-17) 


is  the  time-invariant,  space-variant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz, 
gx(f,h0],hl2)  is  the  scattering  amplitude  function  of  the  discrete  point  scatterer  in  meters, 


”0.1  = 


(3.1-18) 


and 


”l,2  = 


1 1,2 


ll,2 


(3.1-19) 


are  the  dimensionless  unit  vectors  in  the  directions  measured  from  the  source  to  the  discrete  point 
scatterer,  and  from  the  discrete  point  scatterer  to  the  receiver,  respectively,  «(/)  is  the  real, 
frequency-dependent,  attenuation  coefficient  in  nepers  per  meter. 


(3.1-20) 


is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin 
to  appear  at  the  receiver),  and 


ro,i  = 


(3.1-21) 


and 


ri,2 — 


(3.1-22) 


are  the  distances  (ranges)  in  meters  measured  from  the  source  to  the  discrete  point  scatterer,  and 
from  the  discrete  point  scatterer  to  the  receiver,  respectively. 
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The  general  relationship  between  the  acoustic  pressure  pit,  r)  in  pascals  (Pa)  and  the  velocity 
potential  yM(t,  r)  in  squared-meters  per  second  is  given  by 


(3.1-23) 


where  p0(r)  is  the  ambient  (equilibrium)  density  of  the  ocean  in  kilograms  per  cubic  meter,  shown 

as  a  function  of  position  r  =  (x,y,z).  Note  that  lPa  =  lN/m2  where  1 N  =  1  kg -m /sec2. 
Therefore,  upon  substituting  (3.1-16)  into  (3.1-23),  and  since  the  ambient  density  is  constant  in  our 
problem,  the  acoustic  pressure  in  pascals  incident  upon  the  receiver  at  r2  =  (x2,y2,z2)  is  given  by 


pit, r2)  =  -j2nfp0S0HM(f,r2 \r0)exp(+j2nft),  t  >  x. 


(3.1-24) 


where  HM(  f,r2 1  r0 )  is  given  by  (3.1-17).  Note  that  (3.1-24)  can  be  rewritten  as 

p(t,  r2)  =  pf  (r2 )  exp(+j2nft),  t>  r ,  (3.1  -25) 

where  the  spatial-dependent  part  of  the  time-harmonic  acoustic  pressure  in  pascals  is  given  by 

Pf  (r2  )  =  Po  So  H M  (  f  1 1*2  1 1*0  )  •  (3-1  -26) 

Let  us  next  relate  the  scattering  amplitude  function,  the  differential  scattering  cross  section, 
and  target  strength  of  the  discrete  point  scatterer.  The  target  strength  (TS)  is  defined  as  follows 
[13]: 


TS  =  101og10 

Vd{f>K  xAi) 

^ref 

dB  re  AKf, 

(3.1-27) 


where  [12, 13] 


°d(fA b.i Aa)-  lim 

'  '  ru~*° 


r\,2  ^avgs  (  r2  ) 

_  gl(/^0,O«l,2)| 

(4^)2 

(3.1-28) 


is  the  differential  scattering  cross  section  with  units  of  squared  meters,  Iavgj  (rt )  and  Iavg  (r2)  are  the 
time-average,  incident  and  scattered  intensities,  respectively,  with  units  of  watts  per  squared  meter, 
Sif/AiAa)  is  the  scattering  amplitude  function  of  the  discrete  point  scatterer  with  units  of 
meters,  and  Anf  is  a  reference  cross-sectional  area  commonly  chosen  to  be  equal  to  1  m2 . 
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Example  3.1-1  Monostatic  Scattering  Geometry 

For  a  monostatic  (backscatter)  scattering  geometry,  both  the  transmitter  and  receiver  are 
located  at  the  same  position,  that  is. 


Therefore  (see  Fig.  3.1-1), 


and 


r2  =  r0. 


ri,2  ~  r0,l 


(3.1-29) 


(3.1-30) 


«1.2=  “  ”0.1- 


(3.1-31) 


With  the  use  of  (3.1-30)  and  (3.1-31),  the  time-invariant,  space-variant,  complex  frequency 
response  of  the  ocean  given  by  (3.1-17)  reduces  to 


/  ,  ,  .  exp[-2a(/)r01] 

( /, : r2 1 r0 )  =  g , ( /„ h0A, -i n0A ) - ,2  ~ txp(-j2jcff ) , 


(4*ro.i) 


(3.1-32) 


where  fi0  l  is  given  by  (3.1-18), 


is  the  time  delay  in  seconds,  and  r0 1  is  given  by  (3.1-21). 


(3.1-33) 


Example  3.1-2  Pulse  Propagation 

In  this  example  we  will  demonstrate  how  to  use  the  pulse-propagation  coupling  equations 
given  by  (2.2-44)  through  (2.2-47)  for  the  bistatic  scattering  problem  shown  in  Fig.  3.1-1.  We 
begin  by  evaluating  (2.2-47),  which  is  repeated  below  for  convenience: 

y„(?,rRm )  =  X  J_”  X(f)cTi  {f  +  fc)ST  (f  +  fe)Hu(  t, rRn  | /  +  fc,rT^exp(+j2xft)tf, 


(2.2-47) 

where  yM(t,  rR  )  is  the  complex  envelope  of  the  acoustic  signal  (velocity  potential)  incident  upon 
element  n  in  the  receive  array  at  time  t  and  position  rR^  =  {xRii,yRn,zR)  with  units  of  squared- 
meters  per  second. 

The  complex  frequency  response  of  the  ocean  that  describes  the  bistatic  scattering  problem 
shown  in  Fig.  3.1-1  is  given  by  (3.1-17)  and  is  time-invariant  and  space-variant.  Therefore  [see 
(2.1-29)], 

HM(t,r\f,r0)  =  HM(f,r\r0),  (3.1-34) 
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and  as  a  result. 


H„(t,rK]f  +  fr,rr)^Hu[f  +  fc,rK_\rTiy  (3.1-35) 

Substituting  (3.1-17)  into  (3.1-35)  yields 

/i  \  /  Xexp-a(/  +  fe)(rT,  +  r1R  )  r 

Hm[  t,  rRn  |  /  +  fc,rTi  )  =  gi(f  +  fc’KvKn, ) — " — 16  2  r  - —  exp[-y2^(/  +  fc)TTiR 

7],1  l,/fn 


where  [see  (3.1-18)] 


and  [see  (3.1-19)] 


(3.1-36) 


(3.1-37) 


(3.1-38) 


are  the  dimensionless  unit  vectors  in  the  directions  measured  from  element  i  in  the  transmit  array  to 
the  discrete  point  scatterer,  and  from  the  discrete  point  scatterer  to  element  n  in  the  receive  array, 
respectively. 


rT,.l  +  rl.R. 


' 'T„Rn 


(3.1-39) 


is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  acoustic  signal  transmitted  from 
element  i  in  the  transmit  array  to  begin  to  appear  at  element  n  in  the  receive  array)  [see  (3.1-20)], 
and  [see  (3.1-21)] 


and  [see  (3.1-22)] 


(3.1-40) 


ri .*»  “ 


(3.1-41) 


are  the  distances  (ranges)  in  meters  measured  from  element  i  in  the  transmit  array  to  the  discrete 
point  scatterer,  and  from  the  discrete  point  scatterer  to  element  n  in  the  receive  array,  respectively. 
Substituting  (3.1-36)  into  (2.2-47)  yields 


Nr 

yM(t’rRn )  =  X  (/)exp[+j27r/(r-  tTt  R' ) j df  exp(-j2jtfczTj  ^ ),  (3.1-42) 

i=l 


or 
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(3.1-43) 


(t-TTR)exp(-j2nfcTTiiRn 


), 


where 


/  A  A  ,exp \-a(f  +  fc)(rTl  +  rlR) 

X' Rit(/)  =  Z(/)c7;(/  +  /c)57;.(/  +  /c)g1(/  +  /c,n7:,1,«,/eJ -  2  — - 1 

7),1 


(3.1-44) 


and 


Jc'r^n  (0  =  F/1  |x' ,*o  ( / )}  =  J~  X'Ti ^  ( /) exp(+;2^/r)  df .  (3.1  -45) 

Equation  (3.1-43)  indicates  that  the  complex  envelope,  yM(f,rfin),  of  the  acoustic  signal  (velocity 

potential)  incident  upon  element  n  in  the  receive  array  is  equal  to  the  sum  of  time-delayed,  modified 
versions  of  the  complex  envelopes  of  the  transmitted  electrical  signals  from  each  element  in  the 
transmit  array.  Equation  (3.1-44)  shows  how  the  frequency  spectrum  of  the  complex  envelope  of 

the  transmitted  electrical  signal,  X(/),  gets  modified  by  various  other  important  frequency- 
dependent  functions  in  order  to  produce  the  complex  envelope,  yM(t,rR^ ),  of  the  acoustic  signal 

(velocity  potential)  incident  upon  element  n  in  the  receive  array. 

If  we  next  substitute  (3.1-43)  into  (2.2-46),  then  we  obtain 


yM(t>rKm )  =  X  Refe*»  ~ rTt,R)&xp[+j2nfc(t- r^)]}. 

i=l  1 

Since  x(t)  given  by  (2.2-28)  and  x(t)  given  by  (2.2-30)  are  related  by 

x(t)  =  Re{x(t)exp(+j2nfct)}, 

if  we  let 

x'Ti,Rn  (0  =  Re{*r,tf„  (0exp(+j27r/cr)}, 
then  (3.1-46)  can  be  rewritten  as 

Nr 

=  S  XTi>Rn  ^  ~ 

i=  1 

Substituting  (3.1-49)  into  (2.2-45)  yields 

YM{rt,vR)  =  2  X'T,R{r])^{-j27zr]xTiiR), 


(3.1-46) 

(3.1-47) 

(3.1-48) 

(3.1-49) 

(3.1-50) 
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and  upon  substituting  (3.1-50)  into  (2.2-44),  we  obtain 


y(t,  r^)  J*^2^^^(^)exp|+j'27r7](t  n  —  \,2,...,NR,  (3.1-51) 

1=1 


or,  finally, 


where 


and 


Nt 


y(t,TR  )  ^  XTi,R„(t  7 Ti.R „)»  W  1,2,. 


Xri,RSV)  =  XlRn(V)cRn(V)SRn(V) 
x"'R(t)  =  F;1{x"'R'(71)}  =  [jl'tWexpi+jlTTVOdT]. 


(3.1-52) 

(3.1-53) 

(3.1-54) 


Equation  (3.1-52)  indicates  that  the  output  electrical  signal  y(t,rR  )  from  element  n  in  the  receive 

array  at  time  t  and  position  =  (xR  ,yR  ,zR )  with  units  of  volts  is  equal  to  the  sum  of  time- 

delayed,  modified  versions  of  the  transmitted  electrical  signals  from  each  element  in  the  transmit 
array. 


3.2  Discrete  Point  Scatterer  In  Motion 

In  this  section  we  will  analyze  the  bistatic  scattering  problem  shown  in  Fig.  3.2-1.  The 
transmitter  and  receiver  are  not  in  motion.  Only  the  discrete  point  scatterer  is  in  motion.  For 
example  purposes,  all  three  platforms  are  placed  in  deep  water  so  that  the  ocean  medium  can  be 
treated  as  being  unbounded ,  that  is,  there  are  no  boundaries  and,  hence,  no  reflections.  The  speed  of 
sound  and  ambient  density  are  constants.  Although  the  propagation  of  sound  between  the  source 
and  discrete  point  scatterer,  and  between  the  discrete  point  scatterer  and  receiver  can  be  treated  as 
transmission  through  linear,  time- variant,  space-invariant  filters;  the  overall  solution  for  this  bistatic 
scattering  problem  corresponds  to  transmission  through  a  linear,  time-variant,  space-variant  filter. 
The  presence  of  a  discrete  point  scatterer  in  an  unbounded,  homogeneous  fluid  medium  (i.e.,  a  fluid 
medium  with  constant  speed  of  sound  and  ambient  density)  causes  the  medium  to  be  space- variant. 

Let  the  source  distribution  xM(t,  r)  be  a  motionless,  time-harmonic,  point  source  with  units  of 
inverse  seconds,  that  is,  let  [see  (2.1-8)  and  (2.1-1 1)] 

xu(t,r)  =  S08(r  -  r0)txp(+j2nft),  (3.2-1) 

where  S0  is  the  source  strength  in  cubic  meters  per  second  and  the  impulse  function  <5(r  -  r0) ,  with 
units  of  inverse  cubic  meters,  represents  a  point  source  at  r0  =  (jc0,y0,z0).  The  velocity  vector  of  the 
discrete  point  scatterer,  V] ,  is  given  by 


V,  =  ^v 


(3.2-2) 
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where  Vl  is  the  speed  in  meters  per  second  and  nVi  is  the  dimensionless  unit  vector  in  the  direction 
of  V,.  The  velocity  vector  given  by  (3.2-2)  is  constant,  that  is,  both  the  speed  and  direction  are 
constants  -  there  is  no  acceleration.  Motion  begins  at  time  t  =  t0  =  0  seconds. 

Since  the  discrete  point  scatterer  is  now  in  motion,  the  position  vector  from  the  origin  to  the 
discrete  point  scatterer,  denoted  by  ^(t),  is  a  function  of  time  given  by 

&,(t)  =  r,  +  A/Vp  r>0,  (3.2-3) 

where  r,  =  (x,,}^.?,)  is  the  position  vector  from  the  origin  to  the  discrete  point  scatterer  when 
motion  begins  at  time  t  =  t0=  0  seconds  (see  Fig.  3.2-1),  and 


A  t  =  t-t0,  t>t0, 

=  t,  t>  0. 


(3.2-4) 


Note  that  S’,  (0)  =  r, . 

When  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at  some 
time  t'  seconds  where  t'  >  0,  the  position  vector  from  the  origin  to  the  discrete  point  scatterer  is 
given  by  [see  (3.2-3)  and  Fig.  3.2-2] 

r,'  =  S,(r')  =  r,  +  At' V,,  t'>  0,  (3.2-5) 


where 


A  t'=  t'-t0, 
=  t\ 


t'>t0, 
t'>  0. 


(3.2-6) 


The  propagation  of  sound  between  the  source  and  discrete  point  scatterer  can  be  modeled  as 
transmission  through  a  linear,  time-variant,  space-invariant  filter.  Therefore,  the  acoustic  field 
(velocity  potential)  incident  upon  the  discrete  point  scatterer  at  time  t'  and  position  r/=  (x,', y'{, z[ ) 
is  given  by  [see  (2.1-62)] 


yM(t',r;)=  S0HM(t',r;-r0 \f)txp(+j2nft') 

exp(-7^1ri,-rol) 

47r|r1/-r0| 


=  -5fi 


txpi+jlnft'),  t'>  0, 


(3.2-7) 


where 


1/)  =  - 


exp(-yA:|r;-r0|) 
4^|r'- r0| 


(3.2-8) 


By  referring  to  Fig.  3.2-2,  we  can  express  the  position  vector  from  the  point  source  to  the 
discrete  point  scatterer  at  time  t'  as 


ro,i =  ri  “  ro> 


(3.2-9) 
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and  upon  substituting  (3.2-5)  into  (3.2-9),  we  obtain 

ro.i  =  <i«o,i  =  ro,i  +  A''Vr 

where  r0  l  is  given  by  (3.1-4).  Therefore,  (3.2-7)  and  (3.2-8)  can  be  rewritten  as 


and 


yM(t'’ri')=  So  HM{t', r'0)  |  /)  exp(+j2nft' ) 


=  -Sn 


exp(-/Xi 

ti) 

An 

r' 

*0 

i| 

exp(+;'2^/r'),  t'>  0, 


(3.2-10) 


(3.2-11) 


(3.2-12) 


In  order  to  compute  the  acoustic  signal  incident  upon  the  receiver,  we  treat  the  discrete  point 
scatterer  at  time  t>t'  and  position  rf  as  another  motionless,  time-harmonic ,  point  source  with 
units  of  inverse  seconds,  that  is,  let  [see  (3.2-1)  and  Fig.  3.2-2] 

x'M (t, r)  =  S'0S(r  -  if )  exp(+y 2^/0 ,  (3.2-13) 


where  S'0  is  the  source  strength  in  cubic  meters  per  second  and  the  impulse  function  S(r  —  r{), 
with  units  of  inverse  cubic  meters,  represents  a  point  source  at  r'  =  (x[,y[,z[ ) .  The  source  strength 
S'0  is  given  by 


Si  =  S0H„(t',r'0J\f)g,  (f.nij,nU),  (3.2-14) 

where  g1(/,no  I,n1'2)  is  the  scattering  amplitude  function  of  the  discrete  point  scatterer  with  units 

of  meters  as  discussed  in  Section  3.1.  The  unit  of  meters  for  the  scattering  amplitude  function 
represents  an  effective  scattering  length  that  may  be  larger  or  smaller  than  the  actual  length  of  the 
scatterer.  As  we  discussed  in  Section  3.1,  the  scattering  amplitude  function  is  a  complex  function 
(magnitude  and  phase)  and  is,  in  general,  a  function  of  frequency,  the  direction  of  wave  propagation 
from  the  source  to  the  scatterer,  and  the  direction  of  wave  propagation  from  the  scatterer  to  the 
receiver  [12].  Since  the  discrete  point  scatterer  is  now  in  motion  and  the  transmitted  acoustic  field  is 
first  incident  upon  the  discrete  point  scatterer  at  time  t' ,  the  direction  of  wave  propagation  from  the 
source  to  the  scatterer  is  given  by  the  dimensionless  unit  vector  h'0  l  (see  Fig.  3.2-2).  Similarly, 
since  the  discrete  point  scatterer  is  being  treated  as  another  point  source  at  time  t'  and  position  r', 
the  direction  of  wave  propagation  from  the  scatterer  to  the  receiver  is  given  by  the  dimensionless 
unit  vector  h'2  (see  Fig.  3.2-2). 

The  propagation  of  sound  between  the  discrete  point  scatterer  and  receiver  can  also  be 
modeled  as  transmission  through  a  linear,  time-variant,  space-invariant  filter.  The  scattered 
acoustic  field  is  first  incident  upon  the  receiver  at  some  time  t  seconds  where  t>t'>  0.  Therefore, 
the  acoustic  field  (velocity  potential)  incident  upon  the  receiver  at  time  t  and  position  r2  =  (x2,y2,z2) 
due  to  a  point  source  at  r,'  =  {x[,y[,z[ )  is  given  by  [see  (2.1-62)] 
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yM(t,  r2)  =  So  H„(t,  r2  -  r/  |/)exp(+/2*r/f) 


(3.2-15) 


4n\r2-r; 


7] — -  exp(+j2jzft) ,  t>t'>  0, 


where 


f.r  r-|f1_  exp(-jfc|r2-r,'|) 

V’  5  1  W"  4^|r,-r;| 


(3.2-16) 


By  referring  to  Fig.  3.2-2,  we  can  express  the  position  vector  from  the  discrete  point  scatterer 
to  the  receiver  at  time  t'  as 


ri,2  =  *2  ~  ri  * 

and  upon  substituting  (3.2-5)  into  (3.2-17),  we  obtain 

r'  =  r'  n'  =  r  -A t'\ 

*1,2  'l,2Ml,2  *1,2  L*1  vl» 

where  r,  2  is  given  by  (3.1-10).  Therefore,  (3.2-15)  and  (3.2-16)  can  be  rewritten  as 

yM(t,r2)  =  S'0  Hm  (r,  r' 2 1  /)  exp(+j2nft) 
exp(-jk\r'\) 

=  _  So -  |  7  j— ’ -  exp(+j2nft) ,  t>t'>  0, 

4^ki,2 


(3.2-17) 


(3.2-18) 


(3.2-19) 


(3.2-20) 


Let  us  now  begin  the  process  of  obtaining  final  expressions  for  both  the  time-harmonic 
velocity  potential  and  acoustic  pressure  incident  upon  the  receiver.  Substituting  (3.2-14)  into  (3.2- 
19)  yields 


yM(t, r2)  =  f,n'0A, h[2 ) HM (r, r'ia | /) exp(+ j2nft) ,  t  >  t'  >  0, 


or,  equivalently. 


where 


yM(t, r2)  =  S0HM(t, r2  |/,r0)exp(+y2^/O,  t>t'>  0, 


(3.2-21) 


(3.2-22) 
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HM(t,r2  \f,r0)=HM (t\ r' , \f)gx(f,n'0A, h[2 )Hm (t, r'2|/) 

=  giif’K.vKi)- 


=  8 1 


exp 

~jk{ 

1*0,1  +  ri,2  ) 

16jz2 

1*0,1  1*1,2 

exp 

-y*(|r,  -  r0  +  Af' V,|  +  |r2  -  r,  -  At'  V,|) 

167T2 

fj  -  r0  +  At '  V,  r2  -  r,  -  At'  Vx 

At '  is  given  by  (3.2-6),  and 


t  =  t'  + 


or 


(3.2-23) 


(3.2-24) 


(3.2-25) 


Note  that  if  the  bistatic  scattering  problem  shown  in  Fig.  3.2-1  corresponded  to  transmission 
through  a  space-invariant  filter,  then  the  complex  frequency  response  given  by  (3.2-23)  would  be  a 
function  of  the  vector  spatial  difference  r2  -  r0,  which  it  is  not. 

In  order  to  simplify  the  complex  frequency  response  given  by  (3.2-23),  we  need  to  derive 
approximate,  simplified  expressions  for  the  ranges  Ir^J  and  Ir^J.  Let  us  begin  with  |r^  ,| .  Since 


lo,i 


=  ro,.  *  r 


0,1’ 


(3.2-26) 


substituting  (3.2-10)  into  (3.2-26)  yields 


ro,i  =  (rw  +  Ar'V,)*(r0,1  +  Ar'VI). 


(3.2-27) 


Expanding  the  right-hand  side  of  (3.2-27)  and  taking  the  square  root  of  both  sides  of  the  resulting 
equation  yields 


or 


1 

r 

/  A 

2" 

II 

1  +  —  (n0,*V1)Af'+l 

VxAt' 

I 

{  ro.i  J 

(3.2-28) 


|ro,i|  =  ro,f 

,  .  V.At' 

1  +  2— - 

’  .  ,  lVxAt'~\ 

(n01*nv,)  + 

1*0,1 

M  O  »* 

^  r0'1  j 

(3.2-29) 
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Equations  (3.2-28)  and  (3.2-29)  are  exact  expressions  for  the  range  |rofl|.  What  we  need  to 
do  next  is  to  derive  a  valid  binomial  expansion  of  the  square-root  factor  in  (3.2-29)  in  order  to 
simplify  the  expression  for  the  range  Ir^J.  Equation  (3.2-29)  can  be  rewritten  as 


where 


=  r0A^jl  +  b  « 


'0,1 


(3.2-30) 


'0,1 


(«0,l#”v,)  + 


1  VtAt' 

2  ro,i 


(3.2-31) 


Note  that  (3.2-30)  is  valid  only  if  |h|<l,  where  b  is  given  by  (3.2-31).  If  we  impose  the  more 
stringent  criterion  that 


'0,1 


(3.2-32) 


or 


(3.2-33) 


then  we  can  use  only  the  first  two  terms  in  the  binomial  expansion.  For  example,  since  the 
maximum  positive  value  of  the  dot  product  term  in  (3.2-31)  is  +1 ,  and  if  (3.2-33)  is  satisfied,  then 

|i>|<  0.1(1 +  0.025)  (3.2-34) 

<  0.1025  «1. 

Therefore,  if  (3.2-33)  is  satisfied,  then  we  can  use  only  the  first  two  terms  in  the  binomial  expansion 
in  (3.2-30).  Doing  so  yields 


r'  ~r  +VA  t' 

10,1  rQ,\^  y\LyL 


(«o,i  *nVl)+  - 


2  r. 


0,1 


(3.2-35) 


Equation  (3.2-35)  can  be  simplified  further  if  we  can  justify  neglecting  the  second  term  inside 
the  square  brackets.  Consider  the  following  criterion:  the  second  term  inside  the  square  brackets  in 
(3.2-35)  can  be  neglected  if 


^10 


1  VjAf 


2  r, 


0,1 


(3.2-36) 


or 
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cos  a 


>5 


v;a  t' 


'0.1 


(3.2-37) 


where  a  is  the  angle  between  the  two  unit  vectors  n01  and  nVi .  The  following  two  ranges  of  values 
for  the  angle  a  will  satisfy  (3.2-37): 


(3.2-38) 


and 


180°  -cos'1 


V  r<u  ) 


<a<180°. 


(3.2-39) 


For  example,  if  =  0.05  [see  (3.2-33)],  then  0°<a<75.5°  and  104.5°  <a<  180°  will 

satisfy  (3.2-37).  Therefore,  if  (3.2-37)  is  satisfied,  then  (3.2-35)  reduces  to 


r0,l  +  («0,l#Vl)A?/- 


(3.2-40) 


Equation  (3.2-40)  is  our  desired  simplified  expression  for  the  range  .  Let  us  next  derive  a 

similar  expression  for  the  range  |r,'2| . 

Since 


1,2 


=  ri,2#ri,2> 


(3.2-41) 


substituting  (3.2-18)  into  (3.2-41)  yields 


r'u2  =  (r12  -  Ar'Vj)»  (rU2  -  At'  V,). 


(3.2-42) 


Expanding  the  right-hand  side  of  (3.2-42)  and  taking  the  square  root  of  both  sides  of  the  resulting 
equation  yields 


<N 

2" 

|ru|  “  ri,2 

1 - (n12*V.)Ar  + 

l  ru  j 

(3.2-43) 


or 
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1/2 


(3.2-44) 


.  . 

|ru|=ri,2- 

1-2-1 - 

fi,2 

(nl2  •  «v, ) 


i  v;a^ 

^  ri,2 


Equations  (3.2-43)  and  (3.2-44)  are  exact  expressions  for  the  range  |r£2|.  What  we  need  to 
do  next  is  to  derive  a  valid  binomial  expansion  of  the  square-root  factor  in  (3.2-44)  in  order  to 
simplify  the  expression  for  the  range  Ir'J.  Equation  (3.2-44)  can  be  rewritten  as 


1,2 


=  ru2-\l ~b~  q,2 


2  8 


\b\<l. 


(3.2-45) 


where 


b  =  2 


VxAt' 


'1,2 


(nl2*nVi)~ 


1 

2  ri,2 


(3.2-46) 


Note  that  (3.2-45)  is  valid  only  if  |*|<1,  where  b  is  given  by  (3.2-46).  If  we  impose  the  more 
stringent  criterion  that 


2^<0.1. 


(3.2-47) 


'1,2 


or 


YiAL< 

ri,2 


0.05, 


(3.2-48) 


then  we  can  use  only  the  first  two  terms  in  the  binomial  expansion.  For  example,  since  the 
maximum  negative  value  of  the  dot  product  term  in  (3.2-46)  is  -1 ,  and  if  (3.2-48)  is  satisfied,  then 


|*|<  0.1(1  +  0.025) 
<  0.1025  «1. 


(3.2-49) 


Therefore,  if  (3.2-48)  is  satisfied,  then  we  can  use  only  the  first  two  terms  in  the  binomial  expansion 
in  (3.2-45).  Doing  so  yields 


r/  ~  jr  —  1/  At' 
rl,2  M,2  V\LAl 


(  «1,2  •«(',)- 


1  VjAf 


2  r, 


1,2 


(3.2-50) 


Equation  (3.2-50)  can  be  simplified  further  if  we  can  justify  neglecting  the  second  term  inside 
the  square  brackets.  Consider  the  following  criterion:  the  second  term  inside  the  square  brackets  in 
(3.2-50)  can  be  neglected  if 
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(3.2-54) 

For  example,  if  V;At'/r12  =  0.05  [see  (3.2-48)],  then  0°<£<75.5°  and  104.5° <p<  180°  will 
satisfy  (3.2-52).  Therefore,  if  (3.2-52)  is  satisfied,  then  (3.2-50)  reduces  to 

(3.2-55) 

Equation  (3.2-55)  is  our  desired  simplified  expression  for  the  range  |r'2j.  We  are  now  in  a  position 
to  simplify  the  complex  frequency  response  given  by  (3.2-23). 

We  begin  the  simplification  of  the  complex  frequency  response  given  by  (3.2-23)  by 

approximating  the  amplitude  factor  l/(|r'  ,||r'2|).  With  the  use  of  (3.2-40)  and  (3.2-55),  we  can 
express  the  ranges  |r^,|  and  | r'2 [  as  follows: 

1  +  - — cos  a  (3.2-56) 

r0A 

and 


(3.2-57) 
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where  a  is  the  angle  between  the  two  unit  vectors  n01  and  hV[ ,  and  ft  is  the  angle  between  the  two 
unit  vectors  nl  2  and  nV| .  With  the  use  of  (3.2-33)  and  (3.2-48),  (3.2-56)  and  (3.2-57)  reduce  to 


and 


respectively.  Therefore, 


1  1 


/ 

/ 

r0.1 

1 1,2 

(3.2-59) 


(3.2-60) 


Although  these  approximations  of  the  ranges  |rotl|  and  |r'2|  are  acceptable  to  approximate  the 
amplitude  factor,  they  are  not  acceptable  to  approximate  the  phase  factor  expj-j7^|rQ  ,|  +  |rli2|jj  in 
(3.2-23). 

Let  us  approximate  the  phase  factor  exp[— yfc(|ro  x|  +  |r^2|)J  next.  Adding  (3.2-40)  and  (3.2- 
55)  yields 


+  r; 


1,2 


“  rw  +  r1>2  +  [( «0,i  “  «i,2 )  *  vi  ]■ '  > 


(3.2-61) 


where  At'  is  given  by  (3.2-6).  Substituting  (3.2-55)  into  (3.2-25),  and  since  At'- 1'  [see  (3.2-6)], 
we  obtain 


At'* 


(i) 

2 


(3.2-62) 


where 


(3.2-63) 


is  the  dimensionless ,  normalized  (normalized  by  the  constant  speed  of  sound  c)  radial-component 
of  V,  in  the  direction  nl  2.  Substituting  (3.2-62)  into  (3.2-61)  yields 


r\2~ 

|r'  1  +  Ir.'J  «  r0 1  +  r,  2  -  (5  -  l)c 

t-  — 

I  U,  1  j  1  1  V,  1 

C 

(3.2-64) 


or 
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(3.2-65) 


ro.i  +  ru  25  ro,i  +  ru ~(s-l)ct  +  (s-l) r,2. 


where 


1-C 

S=1-C 

(3.2-66) 

is  the  dimensionless ,  time-compression  /  time-stretch  factor,  and 

YL  •  V 
Od)  _  "0,1  M 

(3.2-67) 

P  0,1  ~ 

c 

is  the  dimensionless ,  normalized  (normalized  by  the  constant  speed  of  sound  c)  radial-component 
of  V,  in  the  direction  n0  l.  With  the  use  of  (3.2-65),  the  phase  factor  can  be  expressed  as  follows: 


exp  -;*(|r'  I|  +  |r'2|)  «  exp[+;2^/(5  -  l)r]expj  -jlnf 


(3.2-68) 


or 


where 


exp[-y*(|r;,|  +  |r'2|)]  =  exp[+j2nfs(t-r)]exp(-j2jtft). 


(3.2-69) 


T 


1 

S 


To  +  Cs-i) 


(3.2-70) 


is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin 
to  appear  at  the  receiver),  and 


(3.2-71) 


is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmitter  and  receiver 
are  in  motion. 

As  we  discussed  in  Section  3.1,  in  order  to  model  the  effects  of  frequency-dependent 
attenuation,  simply  replace  the  real  wavenumber  k  given  by  (2.1-23)  with  the  following  complex 
wavenumber  K  [see  (3.1-15)]: 
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K  =  k-ja(f), 


(3.2-72) 


where  <x{f)  is  the  real,  frequency-dependent,  attenuation  coefficient  in  nepers  per  meter.  In 
addition  to  being  real  quantities,  both  k  and  a(  f)  are  positive.  Replacing  the  real  wavenumber  k 

in  the  phase  factor  exp[— | Fq  ,  |  + 1 r'  2 1)  with  the  complex  wavenumber  K  given  by  (3.2-72) 

yields 

exp^dr;,!  +  |r'2|)]  =  exp[-v*(|<1|  +  |r'2|)]exp[-a(/)(|r',|  +  |r'2|)],  (3.2-73) 

where  the  phase  factor  exp[-yl(|ro,i|  +  |ru|)]  is  given  by  (3.2-69).  As  far  as  the  decaying 
exponential  expj^-a(/)dro,I|  +  |r'2|j|  is  concerned,  we  use  the  following  approximation  for  the 
distance  Ir^  J  +  Ir'J  rather  than  (3.2-65): 


(3.2-74) 


so  that 


exp[-a(/)(|r'1|  +  |r'2|)  =  exp[-a(/)(r01  +  r12)]. 


(3.2-75) 


Equation  (3.2-74)  can  be  justified  as  follows.  Since  for  all  real-world  problems  Vl  /c  « 1,  the  result 
is  that  |/J^| « 1  and  « 1  [see  (3.2-67)  and  (3.2-63),  respectively].  Therefore,  5  =  1  [see  (3.2- 

66)]  and  as  a  result,  the  third  term  in  (3.2-64)  can  be  neglected  resulting  in  (3.2-74). 

With  the  use  of  (3.2-22),  (3.2-23),  (3.2-60),  (3.2-73),  (3.2-69),  and  (3.2-75);  we  can 
summarize  our  results  as  follows:  for  the  bistatic  scattering  problem  shown  in  Fig.  3.2-1,  the  time- 
harmonic  velocity  potential  in  squared-meters  per  second  incident  upon  the  receiver  at 
r2  =  (x2,y2,z2),  due  to  a  time-harmonic  point  source  at  r0  =  (x0,y0,z0)  and  a  moving  discrete  point 

scatterer  initially  at  r,  =  (*, ,  yx ,  zx ) ,  is  given  by 


yMit, r2)  =  S0HM(t,r2  |/,r0)exp(+y27r/0,  t  >  x. 


where  S0  is  the  source  strength  in  cubic  meters  per  second, 


(3.2-76) 


,(t, r2  |/,r0)  -  g,(  f,n'0.vKi) ^  +  ~ exP[ '+J27tfs(t  ~  r)]cxp(-j2jtft) 

io/r  r0 ,  r1>2 
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is  the  time-variant,  space-variant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz, 
gi(  /,n'  is  the  scattering  amplitude  function  of  the  discrete  point  scatterer  in  meters,  n0{  and 
n'2  are  the  dimensionless  unit  vectors  in  the  directions  measured  from  the  source  to  the  discrete 
point  scatterer,  and  from  the  discrete  point  scatterer  to  the  receiver,  when  the  transmitted  acoustic 
field  is  first  incident  upon  the  discrete  point  scatterer  at  time  t',  respectively,  a(/)  is  the  real, 
frequency-dependent,  attenuation  coefficient  in  nepers  per  meter, 


(3.2-78) 


is  the  dimensionless,  time-compression  /  time-stretch  factor  [see  (3.2-66)], 


(3.2-79) 


is  the  dimensionless,  normalized  (normalized  by  the  constant  speed  of  sound  c)  radial-component 
of  V,  in  the  direction  h01  [see  (3.2-67)], 


»I.2*V1 


c 


(3.2-80) 


is  the  dimensionless,  normalized  (normalized  by  the  constant  speed  of  sound  c)  radial-component 
of  V,  in  the  direction  hl  2  [see  (3.2-63)], 


T 


1 

S 


T0  +  (S-1) 


(3.2-81) 


is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin 
to  appear  at  the  receiver)  [see  (3.2-70)], 


(3.2-82) 


is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmitter  and  receiver 
are  in  motion  [see  (3.2-71)],  and 


'0,1  ‘ 


(3.2-83) 


and 


rL2  = 


ri 
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(3.2-84) 


are  the  distances  (ranges)  in  meters  measured  from  the  source  to  the  discrete  point  scatterer,  and 
from  the  discrete  point  scatterer  to  the  receiver,  respectively  [see  (3.1-21)  and  (3.1-22)],  when 
motion  begins  at  time  t  =  tQ  =  0  seconds. 

In  order  to  evaluate  the  scattering  amplitude  function  £i(/,«o,i  we  neec*  to  derive 
expressions  for  the  unit  vectors  n'0l  and  h{2 .  In  general, 


y' 

n'  =  0A 
'*0,1  / 

|r0,l| 

(3.2-85) 

and 

_  <2 
**1,2  ~  1  ,  1 » 

|ri,2 1 

(3.2-86) 

where  [see  (3.2-10)] 

ro,i =  *0,1  **0,1  =  ro,i  "**  ^  ^1 

(3.2-87) 

and  [see  (3.2-18)] 

ri,2  =  *"l,2  **1,2  =  ri,2  ~  Vj , 

(3.2-88) 

where 

II 

(3.2-89) 

and 

rf  =  [r/  1 
'1,2  |a1,2|* 

(3.2-90) 

Recall  that  the  parameter 

A /'  is  given  by  [see  (3.2-6)] 

A  t'=t'-t0,  t'>t0, 

=  t',  t'>  0, 

(3.2-91) 

where  t'  is  the  time  instant  when  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete 
point  scatterer  and  /0  =  0  seconds  is  the  time  instant  when  motion  begins.  What  we  need  to  do  next 
is  to  derive  an  expression  for  A t'  to  be  substituted  into  (3.2-87)  and  (3.2-88)  so  that  the  unit 

vectors  given  by  (3.2-85)  and  (3.2-86)  can  be  computed. 

The  time  delay  r  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin  to 
appear  at  the  receiver)  can  be  expressed  as  (see  Fig.  3.2-2) 


(3.2-92) 
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Since 


r'^cAt', 

substituting  (3.2-55)  and  (3.2-93)  into  (3.2-92)  yields 


Ar'~ 

rT 

i-C 

c  _ 

(3.2-93) 


(3.2-94) 


where  t  is  given  by  (3.2-81),  rl  2  is  given  by  (3.2-84),  and  P\1}  is  given  by  (3.2-80).  Equation  (3.2- 

94)  is  the  desired  expression  for  At'  that  will  allow  us  to  compute  the  unit  vectors  according  to 
(3.2-85)  through  (3.2-90).  Note  that  (3.2-94)  can  also  be  obtained  directly  from  (3.2-62)  as 
follows.  If  time  t  in  (3.2-62)  is  set  equal  to  the  time  instant  when  the  scattered  acoustic  field  begins 
to  appear  at  the  receiver,  that  is,  if  t  =  t,  then  (3.2-62)  reduces  to  (3.2-94).  The  parameter  At'  given 
by  (3.2-94)  is  a  constant  whereas  At'  given  by  (3.2-62)  is  a  function  of  time  t  where  t>  x . 

The  general  relationship  between  the  acoustic  pressure  p(t,  r)  in  pascals  (Pa)  and  the  velocity 
potential  yu{t, r)  in  squared-meters  per  second  is  given  by 


(3.2-95) 


where  p0( r)  is  the  ambient  (equilibrium)  density  of  the  ocean  in  kilograms  per  cubic  meter,  shown 

as  a  function  of  position  r  =  (x,y,z).  Note  that  lPa  =  lN/m2  where  1 N  =  1  kg -m /sec2. 
Substituting  (3.2-77)  into  (3.2-76)  yields 


„  \  c  I***  -/  \exP[-a(/X'"o,i  +  'h)]  r  •->  £  U  Nl 

yM(t,r2)  =  S0  gt{f,n0l,nl2  J - L— — - i-txp[+j2nfs(t  -  r)j. 


16tt  r01  Tj  2 


t>r. 


(3.2-96) 

Substituting  (3.2-96)  into  (3.2-95),  and  since  the  unit  vectors  n'0  l  and  h[2  are  constants,  and  the 
ambient  density  is  constant  in  our  problem,  the  acoustic  pressure  in  pascals  incident  upon  the 
receiver  at  r2  =  (x2,y2,z2)  is  given  by 


p(t,r2)  =  -j27tfsp0S0HM(t,r2 \f,r0)exp(+j2nft),  t  >  r, 


(3.2-97) 


where  HM(t, r2  |/,r0)  is  given  by  (3.2-77). 

Let  us  next  relate  the  scattering  amplitude  function,  the  differential  scattering  cross  section, 
and  target  strength  of  the  discrete  point  scatterer.  The  target  strength  (TS)  is  defined  as  follows 
[13]: 
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.2-98) 

where  [12, 13] 

(3.2-99) 

is  the  differential  scattering  cross  section  with  units  of  squared  meters,  Iavg.(r()  and  Iavg  ( r2)  are 
the  time-average,  incident  and  scattered  intensities,  respectively,  with  units  of  watts  per  squared 
meter,  ^1(/,n',,n1'2)  is  the  scattering  amplitude  function  of  the  discrete  point  scatterer  with  units 
of  meters,  and  Ald  is  a  reference  cross-sectional  area  commonly  chosen  to  be  equal  to  1  m2 . 
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(3.2-106) 


is  the  dimensionless,  time-compression  /  time-stretch  factor, 


3(1)  _  »0.1  *  ^1 


fl(1)  =  ■ 

A'o.i 


(3.2-107) 


is  the  dimensionless,  normalized  radial-component  of  V,  in  the  direction  n0 ,, 


T  =  -  T0  +  (s  —  1)- 


(3.2-108) 


is  the  time  delay  in  seconds. 


(3.2-109) 


is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmitter  and  receiver 
are  in  motion,  and  r0 ,  is  given  by  (3.2-83). 

Example  3.2-2  No  Motion  -  Bistatic  And  Monostatic  Scattering  Geometries 

We  begin  this  example  by  considering  the  case  when  the  discrete  point  scatterer  is  not  in 
motion  and  the  scattering  geometry  is  bistatic.  Therefore, 


^  =  0, 


(3.2-110) 


and  as  a  result. 


s  =  l. 


(3.2-111) 


Furthermore  (see  Fig.  3.2-2), 


r0,l  =  r0.1> 


(3.2-112) 


«0,1  —  «0,1  ’ 


(3.2-113) 


(3.2-114) 


n\,2  —  ni,2  • 


(3.2-115) 


With  the  use  of  (3.2-111),  (3.2-113),  and  (3.2-115),  the  time-variant ,  space-variant,  complex 
frequency  response  of  the  ocean  given  by  (3.2-77)  reduces  to  the  exact  time-invariant,  space- 
variant,  complex  frequency  response  of  the  ocean  given  by  (3.1-17),  that  is. 
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H„(t, r2  |/,r0)  =  Hm(  /, r2 1 r0)  =  g,(f,nox,h  1>2) 


exp[-«(/)(r01  +  r12)] 
16^2r01ru 


exp(-j2nfz). 


(3.2-116) 


where  n0 ,  and  n)2  are  given  by  (3.1-18)  and  (3.1-19),  respectively, 

_  r0.1  +  ri.2 

c 


(3.2-117) 


is  the  time  delay  in  seconds,  and  r0 ,  and  rl  2  are  given  by  (3.2-83)  and  (3.2-84),  respectively. 

If  the  discrete  point  scatterer  is  not  in  motion  and  the  scattering  geometry  is  monostatic,  then 
with  the  use  of  (3.2-101)  and  (3.2-102),  the  time-invariant,  space-variant,  complex  frequency 
response  of  the  ocean  given  by  (3.2-1 16)  reduces  to  the  exact  time-invariant,  space- variant,  complex 
frequency  response  of  the  ocean  given  by  (3.1-32),  that  is, 


Hu{t,  r2  |/,r0)  =  Hm{  /,  r2 1  r0)  =  g,(f,n  0A,-n0l) 


exp[-2a(/)r01] 


(4^r0j,) 


exp(-j2nft). 


(3.2-118) 


where  n0 ,  is  given  by  (3.1-18), 


c 

is  the  time  delay  in  seconds,  and  r01  is  given  by  (3.2-83). 


(3.2-119) 


Example  3.2-3  Pulse  Propagation 

In  this  example  we  will  demonstrate  how  to  use  the  pulse-propagation  coupling  equations 
given  by  (2.2-44)  through  (2.2-47)  for  the  bistatic  scattering  problem  shown  in  Fig.  3.2-1.  We 
begin  by  evaluating  (2.2-47),  which  is  repeated  below  for  convenience: 

yM^,rR)  =  ^\[j{f)cTi{f  +  fc)STi{f  +  fc)HM{Kt,vR\f^fc,vT^txp{+j2jzft)df, 

«=i 

(2.2-47) 

where  yM(t, rR  )  is  the  complex  envelope  of  the  acoustic  signal  (velocity  potential)  incident  upon 
element  n  in  the  receive  array  at  time  t  and  position  rR„  =(xR^,yR  ,zRii)  with  units  of  squared- 
meters  per  second. 

The  time-variant,  space-variant,  complex  frequency  response  of  the  ocean  that  describes  the 
bistatic  scattering  problem  shown  in  Fig.  3.2-1  is  given  by  (3.2-77).  Therefore 
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/  |  \  /  Xexp  \-cc(f  +  f)(rT  ,  +  rlR  )  r 

t,rR' \f  +  fc,rTi)~gl(f+ fcJ'Tl.A. )  1  16  2  / - '-^txV[+j2n{f  +  fc)sTi^t 

7},1  IR„ 

xexp[-7'2^(/-f/c)r], 


-r. 


*!■*« 


where 


ST„Rn~' 


1-PZ 


1  -  0(1) 

is  the  dimensionless,  time-compression  /  time-stretch  factor  [see  (3.2-78)], 

pm 


(3.2-120) 


(3.2-121) 


(3.2-122) 


is  the  dimensionless,  normalized  radial-component  of  Vi  in  the  direction  hTj  ,  [see  (3.2-79)  and 
(3.1-37)], 


<«  _  wZl 


pm 


(3.2-123) 


is  the  dimensionless,  normalized  radial-component  of  V,  in  the  direction  hl  R  [see  (3.2-80)  and 
(3.1-38)], 


1 


L 


T0  (STi,R„  1)' 


(3.2-124) 


is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  acoustic  signal  transmitted  from 
element  i  in  the  transmit  array  to  begin  to  appear  at  element  n  in  the  receive  array)  [see  (3.2-81)], 


*o  = 


rT„l  ri,R„ 


(3.2-125) 


is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmit  and  receive 
arrays  are  in  motion  [see  (3.2-82)],  and 


rTi,\—  r7],l 


=  ri-rd 


and 


run.  = 


rUR.  = 


rR  -r, 


(3.2-126) 


(3.2-127) 
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are  the  distances  (ranges)  in  meters  measured  from  element  i  in  the  transmit  array  to  the  discrete 
point  scatterer,  and  from  the  discrete  point  scatterer  to  element  n  in  the  receive  array,  respectively 
[see  (3.2-83)  and  (3.2-84)],  when  motion  begins  at  time  t  =  t0  =  0  seconds. 

In  addition,  in  order  to  evaluate  the  scattering  amplitude  function,  we  need  the  following 
equations: 


and 


nT„l  ~ 


—  *212 


7;.,i 


ni,R~ 


ll  ,R„ 


L1  ,Rn 


(3.2-128) 


(3.2-129) 


are  the  dimensionless  unit  vectors  in  the  directions  measured  from  element  i  in  the  transmit  array  to 
the  discrete  point  scatterer,  and  from  the  discrete  point  scatterer  to  element  n  in  the  receive  array, 
when  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at  time  t  , 
respectively  [see  (3.2-85)  and  (3.2-86)],  where  [see  (3.2-87)] 

r'il  =  r7.,1  +  A?'V1  (3.2-130) 


and  [see  (3.2-88)] 


At'Vj, 


(3.2-131) 


where  rr  ,  =  r,  -  rT  and  r,  R  =rR  -  r,  are  the  position  vectors  measured  from  element  i  in  the 

transmit  array  to  the  discrete  point  scatterer,  and  from  the  discrete  point  scatterer  to  element  n  in  the 
receive  array,  respectively  [see  (3.1-4)  and  (3.1-10)],  when  motion  begins  at  time  t  =  t0  =  0  seconds, 
where  [see  (3.2-94)] 


At'  ~ 


ZT,,R„ 


(3.2-132) 


Substituting  (3.2-120)  into  (2.2-47)  yields 

yM(t,rR" )  ~  ^  J  „^r;.«„(/)exP  +J2jrfsTi.R„(t  ~  tt,,r„ )] #exp[+7'2^/c  sT,,R„(t~  TR.-R»  )]exP 

(3.2-133) 


i=i 


or 


(3.2-134) 


i=l 
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where 


.,  ^xp[-a(/  +  /t)(r,.J  +  ru!,)] 

X{f)cT(f  +  fc)ST.(f  +  fc)gxyf  +  fc,nTVniR^  \(tft2r  r 

3%!  !,/?„ 


(3.2-135) 


and 


it)  =  F/{Z' ,Rn  (/)}  =  [J'Ti'Rn  (f)exp(+j2itft)df .  (3.2-136) 

Equation  (3.2-134)  indicates  that  the  complex  envelope,  yM(t,rRJ,  of  the  acoustic  signal  (velocity 

potential)  incident  upon  element  n  in  the  receive  array  is  equal  to  the  sum  of  time-compressed  / 
time-stretched,  time-delayed,  Doppler  shifted,  modified  versions  of  the  complex  envelopes  of  the 
transmitted  electrical  signals  from  each  element  in  the  transmit  array.  The  amount  of  Doppler  shift 
in  hertz  is  given  by 


^TiyR„  1  X/c- 


(3.2-137) 


Equation  (3.2-135)  shows  how  the  frequency  spectrum  of  the  complex  envelope  of  the  transmitted 
electrical  signal,  X(f),  gets  modified  by  various  other  important  frequency-dependent  functions  in 
order  to  produce  the  complex  envelope,  yM{t,rR„ ),  of  the  acoustic  signal  (velocity  potential)  incident 
upon  element  n  in  the  receive  array. 

Before  continuing  with  the  analysis,  several  comments  concerning  the  dimensionless,  time- 
compression  /  time-stretch  factor  sTj  R  and  (3.2-137)  are  in  order.  There  are  three  ranges  of  values 

for  sT  R  ;  namely,  sT •  R  =  1,  sTjR>  1,  and  sTtR<  1.  When  sTjtR=  1,  this  is  an  indication  that  there  is 

no  time  compression  or  time  stretch,  and  as  a  result,  there  is  no  change  in  signal  bandwidth.  Also, 
when  sTi  R  =1,  ^>TR=0\ that  is,  there  is  no  Doppler  shift  [see  (3.2-137)].  For  example,  if  there  is 

no  motion,  then  the  velocity  vector  of  the  discrete  point  scatterer  is  equal  to  the  zero  vector. 
Therefore,  both  beta  parameters  are  equal  to  zero  and,  as  a  result,  sT  R=  1.  When  sT  <R>  1,  this  is 

an  indication  of  time  compression,  that  is,  the  duration  of  the  received  signal  given  by  (3.2-134)  is 
decreased  (compared  to  the  duration  of  the  transmitted  signal)  resulting  in  an  increase  in  signal 
bandwidth.  Also,  when  sTR  >  1,  <pTR  >  0;  which  indicates  a  positive  Doppler  shift  [see  (3.2-137)]. 

Finally,  when  sT  R  <  1,  this  is  an  indication  of  time  stretch,  that  is,  the  duration  of  the  received  signal 

is  increased  (compared  to  the  duration  of  the  transmitted  signal)  resulting  in  a  decrease  in  signal 
bandwidth.  Also,  when  sT  R  <1,  <pT  R  < 0;  which  indicates  a  negative  Doppler  shift  [see  (3.2- 

137)]. 

If  we  next  substitute  (3.2-134)  into  (2.2-46),  then  we  obtain 


yM(t,rR  )  ~  y,  Rej-^r.,g  (st . r  [*  ^Ti,s,'i')^x.j^+j27cfcsTitRii(t  t,7.  /f  )j|.  (3.2-138) 


Since  x(t)  given  by  (2.2-28)  and  x(t)  given  by  (2.2-30)  are  related  by 


x(t)  =  Re{x(t)exp(+j2nfct)}, 


(3.2-139) 
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if  we  let 


x't„rSO  =  Re{*r,  <RJt)  exp(+j2nfc  f )} , 
then  (3.2-138)  can  be  rewritten  as 

NT 

yM(t’rRn)  ~  S  ■xr,,R.(,sr/I/{„[^  —  tt,,r ,])• 


(3.2-140) 


(3.2-141) 


i=i 


Substituting  (3.2-141)  into  (2.2-45)  yields 

Nr 


yM(V,rR)  *  §  | — - — \^Ti,Rir)lsTi,R„ ) exP(~j^7trlTTi,R„ ) ’  (3.2-142) 


i=l  S 


T„R„ 


and  upon  substituting  (3.2-142)  into  (2.2-44),  we  obtain 

y(r,r^)  =  J  J"  X"  ^(77)exp[+;2^t7(t-  :iR)]dr],  n  =  1,2 ,...,NR,  (3.2-143) 


or,  finally, 


y(t,rR)  ~  XTi,R„(^  ~  ^Tt,Rr  )’  ^  —  .  •  »  Nr  , 


(3.2-144) 


i=i 


where 


X'ti.rS7!)  —  I  |  ^t,,rS^/St!’k«  ^rS^Sr^V) 


Ti,Rn 


(3.2-145) 


and 

jc"  jjf)  =  F;l{x"<Rn(r})}  =  l^X"<Rn(r])Qxp(+j27zrit)dri.  (3.2-146) 

Equation  (3.2-144)  indicates  that  the  output  electrical  signal  y(t,r*n )  from  element  n  in  the  receive 
array  at  time  t  and  position  =  {xR^,yR^,zR)  with  units  of  volts  is  equal  to  the  sum  of  time- 
delayed,  modified  versions  of  the  transmitted  electrical  signals  from  each  element  in  the  transmit 
array. 


3.3  All  Three  Platforms  In  Motion 

In  this  section  we  will  analyze  the  bistatic  scattering  problem  shown  in  Fig.  3.3-1.  All  three 
platforms  -  the  transmitter,  discrete  point  scatterer,  and  receiver  -  are  in  motion.  For  example 
purposes,  they  are  placed  in  deep  water  so  that  the  ocean  medium  can  be  treated  as  being 
unbounded,  that  is,  there  are  no  boundaries  and,  hence,  no  reflections.  The  speed  of  sound  and 


58 


I 

i 

j 


ambient  density  are  constants.  Although  the  propagation  of  sound  between  the  source  and  discrete 
point  scatterer,  and  between  the  discrete  point  scatterer  and  receiver  can  be  treated  as  transmission 
through  linear,  time-variant,  space-invariant  filters;  the  overall  solution  for  this  bistatic  scattering 
problem  corresponds  to  transmission  through  a  linear,  time-variant,  space-variant  filter.  The 
presence  of  a  discrete  point  scatterer  in  an  unbounded,  homogeneous  fluid  medium  (i.e.,  a  fluid 
medium  with  constant  speed  of  sound  and  ambient  density)  causes  the  medium  to  be  space-variant. 

The  velocity  vectors  of  the  transmitter,  V0,  the  discrete  point  scatterer,  V,,  and  the  receiver,  V2, 
are  given  by 

\  =  V0nVo,  (3.3-1) 

V,  =  V;nv  (3.3-2) 

and 

V2=V2nv  (3.3-3) 

where  V0,  Vv  and  V2  are  the  speeds  in  meters  per  second  of  the  transmitter,  the  discrete  point 
scatterer,  and  the  receiver,  respectively,  and  nVo ,  nv, ,  and  hVi  are  the  dimensionless  unit  vectors  in 
the  directions  of  V0,  V,,  and  V2,  respectively.  The  velocity  vectors  given  by  (3.3-1)  through  (3.3-3) 
are  constant ,  that  is,  the  speeds  and  directions  are  constants  -  there  is  no  acceleration.  Motion 
begins  at  time  t  =  t0  =  0  seconds. 

Since  all  three  platforms  are  now  in  motion,  the  position  vectors  from  the  origin  to  the 
transmitter,  discrete  point  scatterer,  and  receiver  -  denoted  by  &0(t),  &t(t),  and  &2(t),  respectively  - 
are  functions  of  time  given  by 


^0(O  =  r0  +  ArVo, 

?>0, 

(3.3-4) 

*,(/)  =  r,  +  At  V„ 

r>0, 

(3.3-5) 

^2(0  =  r2+AfV2, 

t>  0, 

(3.3-6) 

where  r0  =  (x0,y0,z0),  r,  =  O^y,,*,),  and  r2  =  (x2,y2,z2)  are  the  position  vectors  from  the  origin  to 
the  transmitter,  discrete  point  scatterer,  and  receiver,  respectively,  when  motion  begins  at  time 
t  =  t0=  0  seconds  (see  Fig.  3.3-1),  and 


A  t~t  V  o>  (33-1) 

=  t,  t>  0. 

Note  that  ^0(0)  =  r0,  #,(0)  =  rp  and  &2(0)  =  r2 .  Instead  of  trying  to  solve  this  bistatic  scattering 
problem  directly  with  all  three  platforms  in  motion,  we  will  first  create  an  equivalent  problem 
involving  the  transmitter  and  the  discrete  point  scatterer  where  we  can  treat  the  transmitter  (sound 
source)  as  being  motionless.  We  will  then  create  a  second  equivalent  problem  involving  the  discrete 
point  scatterer  and  the  receiver  where  we  can  treat  the  discrete  point  scatterer  (acting  as  a  sound 
source)  as  being  motionless.  This  can  be  accomplished  by  working  with  relative  velocity  vectors. 
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When  motion  begins  at  time  t  =  t0  =  0  seconds,  the  scalar  component  of  V0  in  the  direction 
of  V;  is  given  by  (see  Fig.  3.3-2) 

nVt  •  V0  =  fiVi  •  V0hVo  =  V0(hVi  •  nVo ).  (3.3-8) 

Therefore,  the  velocity  vector  of  the  discrete  point  scatterer  relative  to  the  velocity  vector  of  the 
transmitter  in  the  direction  of  the  velocity  vector  of  the  discrete  point  scatterer  hVi  is  given  by 

(3.3-9) 

By  using  the  relative  velocity  vector  V1>0,  the  transmitter  (sound  source)  can  be  treated  as  being 
motionless .  Therefore,  when  motion  begins  at  time  t  =  t0  =  0  seconds,  the  source  distribution 
xM(t,r)  will  be  treated  as  emotionless,  time-harmonic,  point  source  with  units  of  inverse  seconds, 
that  is,  let  [see  (2.1-8)  and  (2.1-11)] 

xM(t,  r)  =  S0<5(r  -  r0)exp(+j2nft),  (3.3-10) 

where  S0  is  the  source  strength  in  cubic  meters  per  second  and  the  impulse  function  <5(r  -  r0) ,  with 
units  of  inverse  cubic  meters,  represents  a  point  source  at  r0  =  (x0,y0,z0).  In  addition,  since  we  are 
now  working  with  the  relative  velocity  vector  V,  0,  we  need  to  introduce  the  new  position  vector 

gt°’0)(t)  =  r,  +  At  V,  0,  t>  0,  (3.3-11) 

where  At  is  given  by  (3.3-7).  Compare  (3.3-11)  with  (3.3-5).  Note  that  ^’1a,0)(0)  =  r,.  Also  note 
that  if  the  transmitter  is  not  in  motion,  then  V0  =  0,  and  as  a  result,  V,  0  =  V,  [see  (3.3-9)]  and 
^n'0)(t)  =  je,(r)  [see  (3.3-11)  and  (3.3-5)].  And  if  the  discrete  point  scatterer  is  not  in  motion,  then 
Vj  =  0  and  the  relative  velocity  vector  V,0  given  by  (3.3-9)  is  undefined.  In  this  case  we  set 

V,  0  =  0 ,  and  as  a  result,  ^l,0)(t)  =  &,(t)  =  r, . 

When  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at  some 
time  t'  seconds  where  t'  >  0,  the  position  vector  from  the  origin  to  the  discrete  point  scatterer  is 
given  by  [see  (3.3-11)  and  Fig.  3.3-3] 

r1/  =  ^r-0)(t/)  =  rI+At'Vli0,  r'>0,  (3.3-12) 

where 


M'=t'-t0,  t'  >  V 

=  t',  t'  >  0. 


(3.3-13) 


The  propagation  of  sound  between  the  source  and  discrete  point  scatterer  can  be  modeled  as 
transmission  through  a  linear,  time-variant,  space-invariant  filter.  Therefore,  the  acoustic  field 
(velocity  potential)  incident  upon  the  discrete  point  scatterer  at  time  t'  and  position  r(  =  {x[,y[,z[ ) 
is  given  by  [see  (2.1-62)] 
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Figure  3.3-3.  Bistatic  scattering  geometry  when  the  transmitted  acoustic  field  is  first  incident 
upon  the  discrete  point  scatterer  at  time  t'  seconds  and  when  the  scattered  acoustic  field  is  first 
incident  upon  the  receiver  at  time  t  seconds  where  t>t'>t0  (t0  =  0).  Point  0,  /J,(r0),  is  the 
transmitter;  point  1,  is  the  discrete  point  scatterer;  and  point  2,  P2(r2),  is  the  receiver.  All 

three  platforms  are  in  motion. 


yM(t',r()  =  S0HM{t',r(-r0 \f)exp(+j2nft') 


(3.3-14) 


where 


=  -S0 — j — 7 - -exp(+y'2;r/r')>  *'  >  0, 

An  r,  -r0 


„  r|fl  exp(-,/*|r,'-r1,|) 

M‘-r<  r«l  f)-  4?dr,'-r,,l  ' 


(3.3-15) 


By  referring  to  Fig.  3.3-3,  we  can  express  the  position  vector  from  the  point  source  to  the 
discrete  point  scatterer  at  time  t'  as 


ro,i  =  ri/-ro’ 


and  upon  substituting  (3.3-12)  into  (3.3-16),  we  obtain 

r0,, =  ro,i «o.i =  ro,i  ^  V1>0, 

where  r0 ,  is  given  by  (3.1-4).  Therefore,  (3.3-14)  and  (3.3-15)  can  be  rewritten  as 
yM(t',<)  =  S0HM(t',r'A\f)txp(+j2nft') 


=  -Sn 


3(~7^|ro,i|) 


I  ’ '  exp(+ j2nft'),  t' >  0, 


(3.3-16) 


(3.3-17) 


(3.3-18) 


exp(~7^|rQ1|j 

4^|ro.i| 


(3.3-19) 


Let  us  now  create  a  similar  equivalent  problem  involving  the  discrete  point  scatterer  and  the 
receiver  where  we  can  treat  the  discrete  point  scatterer  (acting  as  a  sound  source)  as  being 
motionless.  When  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at 
time  t'  >  0  seconds,  the  position  vector  from  the  origin  to  the  receiver  is  given  by  [see  (3.3-6)  and 
Fig.  3.3-3] 


r'  =  Z2(t')  =  r2+At'Y2,  t' >  0, 


(3.3-20) 


where  At'  is  given  by  (3.3-13).  Equation  (3.3-20)  indicates  that  after  At'  seconds,  the  receiver  - 
independent  of  the  transmitter  and  the  discrete  point  scatterer  -  travels  an  additional  distance  of 
A t'V2  meters  in  the  direction  hVi  (see  Fig.  3.3-3).  Also  at  time  t',  the  scalar  component  of  V,  in 

the  direction  of  V2  is  given  by  (see  Fig.  3.3-4) 


ftu  •  Vj  =  nv  •  Vx hv  =  Vx(hv  *nV]). 


(3.3-21) 
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Therefore,  the  velocity  vector  of  the  receiver  relative  to  the  velocity  vector  of  the  discrete  point 
scatterer  in  the  direction  of  the  velocity  vector  of  the  receiver  nVj  is  given  by 


V, ,  =  V2 -(nVi  .  V,)^  =  [v,  -  (((«„,  •  \  )]*„, 


(3.3-22) 


By  using  the  relative  velocity  vector  V2I,  the  discrete  point  scatterer  (acting  as  a  sound  source)  can 
be  treated  as  being  motionless.  Therefore,  in  order  to  compute  the  acoustic  signal  at  the  receiver,  we 
treat  the  discrete  point  scatterer  at  time  t  >  t'  and  position  r'  as  another  motionless,  time-harmonic, 
point  source  with  units  of  inverse  seconds,  that  is,  let  [see  (3.3-10)  and  Fig.  3.3-3] 

x'M(t,r)  =  Sq8(t  -  r')exp(+j2nft),  (3.3-23) 

where  is  the  source  strength  in  cubic  meters  per  second  and  the  impulse  function  <5(r  —  r/), 
with  units  of  inverse  cubic  meters,  represents  a  point  source  at  r(  =  (x[,  y{,z{),  where  is  given  by 
(3.3-12).  The  source  strength  Sq  will  be  given  later. 

Since  we  are  now  working  with  the  relative  velocity  vector  V2 ,,  we  need  to  introduce  the  new 
position  vector 


^2,I)(r)  =  r2  +  At"\u,  t>t'>  0, 


(3.3-24) 


where 


A t"=t-t',  t>t'>  0,  (3.3-25) 

and  r2  is  given  by  (3.3-20).  Compare  (3.3-24)  with  (3.3-6).  Note  that  J£22,1)(0  =  r2 .  Also  note  that 
if  the  discrete  point  scatterer  is  not  in  motion,  then  V,  =  0,  and  as  a  result,  V2 ,  =  V2  [see  (3.3-22)] 
and  jef,l)(r)  =  752{t)  for  t>t'>0  [see  (3.3-24)  and  (3.3-6)].  And  if  the  receiver  is  not  in  motion, 
then  V2  =  0  and  the  relative  velocity  vector  V21  given  by  (3.3-22)  is  undefined.  In  this  case  we  set 

V2 ,  =  0,  and  as  a  result,  ^,22,1)(r)  =  ^2(0  =  r2  for  t  >  t'  >  0. 

When  the  scattered  acoustic  field  is  first  incident  upon  the  receiver  at  some  time  t  seconds 
where  t  >  t'  >  0,  the  position  vector  from  the  origin  to  the  receiver  is  given  by  [see  (3.3-24)  and  Fig. 
3.3-3] 


r2"  =  2?<2’1)(0  =  r'  +  Ar"V2>1,  t>t'>0,  (3.3-26) 

where  At"  is  given  by  (3.3-25).  Equation  (3.3-26)  indicates  that  after  At"  seconds,  the  receiver 
travels  an  additional  distance  of  A t"V2  l  meters  in  the  direction  (see  Fig.  3.3-3).  By  referring  to 
Fig.  3.3-3,  it  can  be  seen  that 

=  (3-3-27) 

where  r"2  is  the  position  vector  between  the  discrete  point  scatterer  and  the  receiver  at  time 
t>t'>  0.  Substituting  (3.3-12)  and  (3.3-26)  into  (3.3-27)  yields 
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(3.3-28) 


<2  =  =  ru  +  A t'(\2  - \0)+A t"\ 


2,1’ 


where  r12  is  given  by  (3.1-10).  Equation  (3.3-28)  can  also  be  obtained  directly  from  Fig.  3.3-3  - 
without  tbe  use  of  (3.3-27)  -  simply  by  using  vector  addition.  Doing  so  yields 


r"  =  r"  n"  =  -  At'  Vw  +  ru  +  Ar' V2  +  A /"  V2>1 
=  rlj2  +  At\  V2  -  V,  0 ) + At"  V2il . 


(3.3-29) 


The  propagation  of  sound  between  the  discrete  point  scatterer  and  receiver  can  also  be 
modeled  as  transmission  through  a  linear,  time-variant,  space-invariant  filter.  Therefore,  the 
acoustic  field  (velocity  potential)  incident  upon  the  receiver  at  time  t  and  position  r2"  =  (x",  y",z") 
due  to  a  point  source  at  r'  =  (x'vy[,z[ )  is  given  by  [see  (2.1-62)] 


yM(t,r")=  S'0HM{t,r"-v'  \f)exp{+j2nft) 

=  ~so — rn; — ^-^exp(+;2^/o,  t>t  > o, 

4^|r2  ~ri  | 


where 


/  .  \  exp(— yA:  I  r2"  —  r'  I ) 

With  the  use  of  (3.3-27),  (3.3-30)  and  (3.3-31)  can  be  rewritten  as 
yM(t,r")  =  So  hm  (t>  Ki  |  /)  exp(+y  2  Tift) 

- exp(+ j2n ft),  t>t'>  0, 


exp(-yl|i 

<2  ) 

Atc 

) 

and 


exp(-;fc|r^|) 

An\<i\ 


The  source  strength  S0'  is  given  by  [see  (3.3-18)] 


(3.3-30) 


(3.3-31) 


(3.3-32) 


(3.3-33) 


(3.3-34) 


where  g{(f,n'0 1,nI,/2)  is  the  scattering  amplitude  function  of  the  discrete  point  scatterer  with  units 
of  meters  as  discussed  in  Section  3.1.  The  unit  of  meters  for  the  scattering  amplitude  function 
represents  an  effective  scattering  length  that  may  be  larger  or  smaller  than  the  actual  length  of  the 
scatterer.  As  we  discussed  in  Section  3.1,  the  scattering  amplitude  function  is  a  complex  function 
(magnitude  and  phase)  and  is,  in  general,  a  function  of  frequency,  the  direction  of  wave  propagation 
from  the  source  to  the  scatterer,  and  the  direction  of  wave  propagation  from  the  scatterer  to  the 
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receiver  [12].  Since  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer  at 
time  t' ,  the  direction  of  wave  propagation  from  the  source  to  the  scatterer  is  given  by  the 
dimensionless  unit  vector  h'0  l  (see  Fig.  3.3-3).  Similarly,  since  the  discrete  point  scatterer  is  being 
treated  as  another  point  source  at  time  t'  and  position  r,',  and  the  scattered  acoustic  field  is  first 
incident  upon  the  receiver  at  time  t  >  t'  >  0,  the  direction  of  wave  propagation  from  the  scatterer  to 
the  receiver  is  given  by  the  dimensionless  unit  vector  h"2  (see  Fig.  3.3-3). 

Let  us  now  begin  the  process  of  obtaining  final  expressions  for  both  the  time-harmonic 
velocity  potential  and  acoustic  pressure  incident  upon  the  receiver.  Substituting  (3.3-34)  into  (3.3- 
32)  yields 

y„(t.  If)  =  S„tf„(<X,|/)«,(  |/)exp(+./2tf/t),  t  >  t’  >  0, 


or,  equivalently, 


yM(t, r")  =  S0HM(t,r"  \f,r0)exp(+j2nft),  t  >  t'  >  0, 


(3.3-35) 

(3.3-36) 


where 


exp 

-jk( 

ro,. 

+  ] 

Cl) 

16tt2 

ro,i 

K: 

J 

exp 

~jk( 

|ri  -  ro  +  At'Vjo 

+ 

r2 -  ri  +  At'(  V2  -  V,  0) + At"  V2  j|)] 

16^2 

ri_ro  +  At'V^o 

r2 

-  r,  +  At'(\2  -  V,  0) + A?" V2 , 

(3.3-37) 


At'  is  given  by  (3.3-13),  Vu0  is  given  by  (3.3-9),  At"  is  given  by  (3.3-25),  V24  is  given  by  (3.3- 
22),  and 


t  =  t'  + 


(3.3-38) 


or 


(3.3-39) 


Note  that  if  the  bistatic  scattering  problem  shown  in  Fig.  3.3-1  corresponded  to  transmission 
through  a  space-invariant  filter,  then  the  complex  frequency  response  given  by  (3.3-37)  would  be  a 
function  of  the  vector  spatial  difference  r2"-r0,  which  it  is  not  [r"  is  given  by  (3.3-26)]. 
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In  order  to  simplify  the  complex  frequency  response  given  by  (3.3-37),  we  need  to  derive 
approximate,  simplified  expressions  for  the  ranges  Ir^J  and  I  r^J.  Let  us  begin  with  r^J.  Since 


lo,i 


=  ro,.*ro.P 


(3.3-40) 


substituting  (3.3-17)  into  (3.3-40)  yields 

|ro,i|2  =  (1*0,1  +  At' V10)  •  (r0,i  +  At'Vlt0). 


(3.3-41) 


Expanding  the  right-hand  side  of  (3.3-41)  and  taking  the  square  root  of  both  sides  of  the  resulting 
equation  yields 


ro,i  —  ro,i 


1  +— («o.,*V1,0)Af/  + 


'0,1 


KqA^ 
V  ro,i  ) 


2-|l/2 


(3.3-42) 


or 


L'  I  _  -  , 

L2^Ar,[ 

|r0,l|  ~  r0,l 

1*0.1 

(n0,i*nVl)  + 


1  KqA t' 


2  r, 


0,1 


1/2 


(3.3-43) 


Equations  (3.3-42)  and  (3.3-43)  are  exact  expressions  for  the  range  |ro  ,|.  What  we  need  to 
do  next  is  to  derive  a  valid  binomial  expansion  of  the  square-root  factor  in  (3.3-43)  in  order  to 
simplify  the  expression  for  the  range  1^,1.  Equation  (3.3-43)  can  be  rewritten  as 


ro,i  |  =  rMVI+S  “  ro.i 


( i  b  b2 

1  + - +  • 

^28 


\b\<l. 


(3.3-44) 


where 


b  =  2 


KoAt' 


'0,1 


(n0,,*nVi)  + 


1  KqA^ 

2  ^o.i 


(3.3-45) 


Note  that  (3.3-44)  is  valid  only  if  \b\<l,  where  b  is  given  by  (3.3-45).  If  we  impose  the  more 
stringent  criterion  that 


2^  <0.1, 


'0,1 


(3.3-46) 


or 
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or 

(3.3-51) 

where  a  is  the  angle  between  the  two  unit  vectors  n0  l  and  nVi .  The  following  two  ranges  of  values 
for  the  angle  a  will  satisfy  (3.3-51): 

(3.3-52) 

and 

(3.3-53) 

For  example,  if  Vu0At'/rM  =  0.05  [see  (3.3-47)],  then  0°<a<75.5°  and  104.5° < a<  180°  will 
satisfy  (3.3-51).  Therefore,  if  (3.3-51)  is  satisfied,  then  (3.3-49)  reduces  to 
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(3.3-54) 


~  r0,t  +  («0.i #  V,0)A/  • 


Equation  (3.3-54)  is  our  desired  simplified  expression  for  the  range  ,| .  Let  us  next  derive  a 

similar  expression  for  the  range  |r"2|. 

Since 


ll,2 


_  r"  «  r" 

1 1,2  *1.2 » 


(3.3-55) 


substituting  (3.3-29)  into  (3.3-55)  yields 

kur  =  [ru+A''(V2-V1.o)+A("V2,].[ru+A('(V2-V,,„)+A»"V2.l].  (3.3-56) 


Expanding  the  right-hand  side  of  (3.3-56)  and  taking  the  square  root  of  both  sides  of  the  resulting 
equation  yields 


rf/  —  r 
* 1,2  'l,2 


1  +  —  [n,.;*(V2-V,t)]A('+ 


'  1,2 


lVj  V,'“IA'  '  +— (nl2» V21)Ar"+ 


+2^I[V2,1.(V2-Vi,<1)]|.  , 


V  'U2 

1/2 


'  1,2 


rvM^2 

\  rW  J 


'1,2 


(3.3-57) 


(3.3-58) 

Equations  (3.3-57)  and  (3.3-58)  are  exact  expressions  for  the  range  |r"J.  What  we  need  to 
do  next  is  to  derive  a  valid  binomial  expansion  of  the  square-root  factor  in  (3.3-57)  in  order  to 
simplify  the  expression  for  the  range  |r"2| .  We  begin  the  simplification  process  by  neglecting  those 

terms  in  (3.3-57)  involving  A t'  and  A t"  squared  and  the  cross  product  term  involving  At' At". 
Doing  so  yields 
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1,2 


=  r,2A/T+fc«r,i2 


(  b  b2  ^ 

1  + - +  •■ 

v  2  8 


,  \b\<  1, 


where 


(3.3-59) 


b  =  —  { [n,.2  •  ( V,  -  V10)]&/'+  (nu  •  V2J)Ar"  } . 

ri,2 


(3.3-60) 


Note  that  the  binomial  expansion  in  (3.3-59)  is  valid  only  if  |&|  <  1,  where  b  is  given  by  (3.3-60). 
Assuming  that  \b\  <  1  and  if  we  use  only  the  first  two  terms  in  the  binomial  expansion  in  (3.3-59), 
then 


1,2 


»  '•u  +  [«,.2*(V2  -  V1c,)]A('-K»,,2» 


(3.3-61) 


r"2  .  We  are  now  in  a  position 


Equation  (3.3-61)  is  our  desired  simplified  expression  for  the  range 

to  simplify  the  complex  frequency  response  given  by  (3.3-37). 

We  begin  the  simplification  of  the  complex  frequency  response  given  by  (3.3-37)  by 

approximating  the  amplitude  factor  1^/^ | i*q  1 1 j r^2 1 ^ .  With  the  use  of  (3.3-54)  and  (3.3-61),  we  can 
express  the  ranges  1^,1  and  |r"2|  as  follows: 


r0,l  *  r0,l 


i+s**: 


cos  a 


'0,1 


(3.3-62) 


where  a  is  the  angle  between  the  two  unit  vectors  n0 ,  and  nVi ,  and 

')  ,  K2*(V2~vi.o)|a<'  ,  (*,.,♦  V2,,)At' 


ri,2  ~ri,2 


1,2 


'1,2 


(3.3-63) 


With  the  use  of  (3.3-47)  in  (3.3-62),  and  upon  neglecting  the  second  and  third  terms  inside  the 
square  brackets  in  (3.3-63),  (3.3-62)  and  (3.3-63)  reduce  to 


and 


respectively.  Therefore, 


ro,i  *  ro.i 


ri,2  ~  ri,2  ’ 


(3.3-64) 


(3.3-65) 


r0,l  ri,2  r0,iri,2 


(3.3-66) 
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Although  these  approximations  of  the  ranges  j  ,  |  and  |r"2|  are  acceptable  to  approximate  the 
amplitude  factor,  they  are  not  acceptable  to  approximate  the  phase  factor  exp[^— |  ,  |  + 1  r"2 1 )]  in 


(3.3-37). 


Let  us  approximate  the  phase  factor  expj-y'^ | r', |  +  |r"2 1 jj  next.  Adding  (3.3-54)  and  (3.3- 
61)  yields 

|r0.i|  +  |ri"|  *  h.i  +  r, . 2  +  { [(«0,i  -  «1.2)  *  Vi.o]  +  («1,2  •  V2)}Ar'+ (nl  2  •  V2J)Ar",  (3.3-67) 


where  At'  is  given  by  (3.3-13)  and  At"  is  given  by  (3.3-25).  Substituting  (3.3-61)  into  (3.3-39), 
and  since  A t'=  t'  [see  (3.3-13)],  we  obtain 


c  c  c 

The  parameter  At" ,  given  by  (3.3-25),  can  be  rewritten  as  follows: 


(3.3-68) 


A r-t-f 

=  t-t0-(t'-t0),  t>t'>  t0. 


(3.3-69) 


or 


A t"=  t-t0-  At',  t>t'>  t0, 

=  t- At',  t>t'>  0, 

where  At'  is  given  by  (3.3-13).  Substituting  (3.3-70)  into  (3.3-68)  yields 


At'  ~ 


r+k22,-/sg,>f1 


where 


(3.3-70) 


(3.3-71) 


(3.3-72) 


is  the  dimensionless ,  normalized  (normalized  by  the  constant  speed  of  sound  c )  radial-component 
of  V2  in  the  direction  hl  2 , 


(3.3-73) 
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is  the  dimensionless,  normalized  (normalized  by  the  constant  speed  of  sound  c )  radial-component 
of  V,  0  in  the  direction  nl  2,  and 


P 


(2,1)  _ 
1,2  — 


*  V2,l 


c 


(3.3-74) 


is  the  dimensionless,  normalized  (normalized  by  the  constant  speed  of  sound  c )  radial-component 
of  V2  j  in  the  direction  n,  2.  Substituting  (3.3-70)  and  (3.3-71)  into  (3.3-67)  yields 


5-1 +02° 

r0,l"*'ri,2  (S  l)c^+  j_^g(2,l)  ri,2’ 


(3.3-75) 


where 


[1-C’l 

[l-ffi1’] 

■“[1-W1+ 

is  the  dimensionless,  time-compression  /  time-stretch  factor,  and 


(3.3-76) 


(3.3-77) 


is  the  dimensionless,  normalized  (normalized  by  the  constant  speed  of  sound  c)  radial-component 
of  V,  0  in  the  direction  n01.  With  the  use  of  (3.3-75),  the  phase  factor  can  be  expressed  as  follows: 


exp[-yl(|r'  ,|  +  |r"|)]  «  &xp[+j2nf(s  -  l)/]exp 


-j2nf 


To  +  ' 


■l  +  Pt 


(2,1) 


f  1,2 


(2,1) 

,2 


(3.3-78) 


or 


where 


exp  -y£(|ro,I|  +  |r"|)]  «  exp[+j‘2^/j(t  -  x)]exp(-j2nft). 


(3.3-79) 


1 

■S  1  +  P\,2  }  ri,2 

T  —  — 
5 

0  1-^2° 

(3.3-80) 
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is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin 
to  appear  at  the  receiver),  and 


(3.3-81) 


is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmitter  and  receiver 
are  in  motion. 

As  we  discussed  in  Section  3.1,  in  order  to  model  the  effects  of  frequency-dependent 
attenuation,  simply  replace  the  real  wavenumber  k  given  by  (2.1-23)  with  the  following  complex 
wavenumber  K  [see  (3.1-15)]: 


K  =  k-ja(f), 


(3.3-82) 


where  a(f)  is  the  real,  frequency-dependent,  attenuation  coefficient  in  nepers  per  meter.  In 
addition  to  being  real  quantities,  both  k  and  a(  f)  are  positive.  Replacing  the  real  wavenumber  k 

in  the  phase  factor  expJ-_/&| |  + 1 r"2 1 jj  with  the  complex  wavenumber  K  given  by  (3.3-82) 
yields 

exp[-jX(|r^|  +  |r"|)]  =  exp[-yk(|r;  ,|  +  |r"2|)]exp[-a(/)(|r'tl|  +  |r"|)],  (3.3-83) 

where  the  phase  factor  exp|-y^ |  r^,  |  + 1 ru |)]  is  giyen  by  (3.3-79).  As  far  as  the  decaying 
exponential  exp|  - a( / ) ^ |  r^, |  + 1 r"2 1 j J  is  concerned,  we  use  the  following  approximation  for  the 
distance  |ro1|  +  |r1,'2|  rather  than  (3.3-75): 

|ro,i| + 1  *1.2 1 25  ro,i +  ru  (3.3-84) 

so  that 


exp  -a(/)(|r'  ,|  +  |r1';2|)]  *  exp[-a(/)(r0I  +  r12)]. 


(3.3-85) 


Equation  (3.3-84)  can  be  justified  as  follows.  Since  for  all  real-world  problems  Vl  0/c « 1, 
V21/c  « 1,  and  V2/c  « 1,  the  result  is  that  |)3o’10)|«l,  |j8^21)|«1,  |)8(1|20)|«1,  and  |^2  «1 
[see  (3.3-77),  (3.3-74),  (3.3-73),  and  (3.3-72),  respectively].  Therefore,  s  ~  1  [see  (3.3-76)]  and 
since  « 1,  the  third  and  fourth  terms  in  (3.3-75)  can  be  neglected  resulting  in  (3.3-84). 

With  the  use  of  (3.3-36),  (3.3-37),  (3.3-66),  (3.3-83),  (3.3-79),  and  (3.3-85);  we  can 
summarize  our  results  as  follows:  for  the  bistatic  scattering  problem  shown  in  Fig.  3.3-1,  the  time- 
harmonic  velocity  potential  in  squared-meters  per  second  incident  upon  the  receiver  at 
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r"  -  (x'fy'fz'j),  due  to  a  moving  time-harmonic  point  source  initially  at  r0  =(x0,y0,z0)  and  a 
moving  discrete  point  scatterer  initially  at  r,  =  (xx ,  y, ,  z, ) ,  is  given  by 


yM(t,x")  =  S0HM(t,r"  |  /,  r0 )  exp(+j2jvft),  t  >  t, 


(3.3-86) 


where  S0  is  the  source  strength  in  cubic  meters  per  second, 


(3.3-87) 


is  the  time-variant,  space-variant,  complex  frequency  response  of  the  ocean  at  frequency  /  hertz, 
8i(  f,KvK'i) is  die  scattering  amplitude  junction  of  the  discrete  point  scatterer  in  meters,  n0A  is 
the  dimensionless  unit  vector  in  the  direction  measured  from  the  source  to  the  discrete  point 
scatterer  when  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point  scatterer,  nl  2  is 
the  dimensionless  unit  vector  in  the  direction  measured  from  the  discrete  point  scatterer  to  the 
receiver  when  the  scattered  acoustic  field  is  first  incident  upon  the  receiver,  oc(  f )  is  the  real, 
frequency-dependent,  attenuation  coefficient  in  nepers  per  meter, 


[1-/C] 

[i-jsS"] 

(3.3-88) 


is  the  dimensionless,  time-compression  /  time-stretch  factor  [see  (3.3-76)], 


P 


(1.0)  _ 
0,1  — 


"o.l  *  Vl.O 


c 


(3.3-89) 


is  the  dimensionless,  normalized  (normalized  by  the  constant  speed  of  sound  c )  radial-component 
of  V,  0  in  the  direction  h0 ,  [see  (3.3-77)], 


P 


(2,1)  _ 
1,2  — 


^1,2  *  ^2.1 


c 


(3.3-90) 


is  the  dimensionless,  normalized  (normalized  by  the  constant  speed  of  sound  c)  radial-component 
of  V2 ,  in  the  direction  hl  2  [see  (3.3-74)], 


”l,2*Vl,0 

c 


(3.3-91) 
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is  the  dimensionless ,  normalized  (normalized  by  the  constant  speed  of  sound  c)  radial-component 
of  V,  0  in  the  direction  nl2  [see  (3.3-73)], 


(3.3-92) 


is  the  dimensionless,  normalized  (normalized  by  the  constant  speed  of  sound  c)  radial-component 


of  V2  in  the  direction  hx 

2  [see  (3.3-72)], 

Vu,  =  Vt -(«,,•  %)»,,= 

'k-wM 

nvx 

(3.3-93) 

is  the  velocity  vector  of  the  discrete  point  scatterer  relative  to  the  velocity  vector  of  the  transmitter  in 
the  direction  of  the  velocity  vector  of  the  discrete  point  scatterer  hVi  [see  (3.3-9)], 

v2J  =  v2-(Vv.»iv,= 

^2  V[(>V2  #  ^  )] 

K 

(3.3-94) 

is  the  velocity  vector  of  the  receiver  relative  to  the  velocity  vector  of  the  discrete  point  scatterer  in  the 
direction  of  the  velocity  vector  of  the  receiver  hVi  [see  (3.3-22)], 


1 

_  J-l  +  fta0  ru 

l  — 

s 

r  c  \ 

(3.3-95) 


is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin 
to  appear  at  the  receiver)  [see  (3.3-80)], 


(3.3-96) 


is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmitter  and  receiver 
are  in  motion  [see  (3.3-81)],  and 


^o.i  =  |ro,i|  =  |ri-ro|  (3.3-97) 

and 

ri,2  =  |ri,2|  =  Ir2  ~  ri|  (3.3-98) 

are  the  distances  (ranges)  in  meters  measured  from  the  source  to  the  discrete  point  scatterer,  and 
from  the  discrete  point  scatterer  to  the  receiver,  respectively  [see  (3.1-21)  and  (3.1-22)],  when 
motion  begins  at  time  t  =  t0=  0  seconds.  Recall  that  if  the  discrete  point  scatterer  is  not  in  motion. 
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then  V,  =  0  and  the  relative  velocity  vector  V10  given  by  (3.3-93)  is  undefined.  In  this  case  we  set 
V,  0  =  0 .  Similarly,  if  the  receiver  is  not  in  motion,  then  V2  =  0  and  the  relative  velocity  vector  V2 , 
given  by  (3.3-94)  is  undefined.  In  this  case  we  set  V2 ,  =  0.  The  position  vector  from  the  origin  to 
the  receiver  when  the  scattered  acoustic  field  is  incident  upon  the  receiver  is  given  by  [see  (3.3-26), 
(3.3-20),  and  (3.3-70)] 


r2"  =  r2  +  At'V2  +  (?-Ar')V2>I,  t>  r, 


(3.3-99) 


where  At'  is  given  by  (3.3-71)  and  the  relative  velocity  vector  V2 ,  is  given  by  (3.3-94). 

In  order  to  evaluate  the  scattering  amplitude  function  gx[f,n we  need  to  derive 
expressions  for  the  unit  vectors  n'0  i  and  h"2 .  In  general, 


ro,i 

(3.3-100) 

™0,1  “ 

<1 

and 

C* 

II 

<2 

Ir"  r 

(3.3-101) 

where  [see  (3.3-17)] 

*0.1  =  <i"o,l  = 

roa  +  At'V.,0 

(3.3-102) 

and  [see  (3.3-28)] 

r"2  =  r~  />"  =  ru + At'(V2  -  V,0)+ At"  V2J, 

(3.3-103) 

where 

<1  =  | 

t'  1 

10,l| 

(3.3-104) 

and 

<2  = 

r"  1 

*  1,2  j  * 

(3.3-105) 

Recall  that  the  parameters  At'  and  At"  are  given  by  [see  (3.3-13)  and  (3.3-70)] 

At'=t'-t0,  t'>t0, 

=  t',  t'>  0, 


(3.3-106) 


and 
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A  t"=  t-t0-  A  t', 
=  t  -  At', 


(3.3-107) 


t  >  t'  >  t0, 
t>t'>  0, 


where  t'  is  the  time  instant  when  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete 
point  scatterer  and  t0  =  0  is  the  time  instant  when  motion  begins.  If  time  t  in  (3.3-107)  is  set  equal 
to  the  time  instant  when  the  scattered  acoustic  field  begins  to  appear  at  the  receiver,  that  is,  if  t  —  x , 
then  (3.3-107) reduces  to 

At"=x- At' .  (3.3-108) 


What  we  need  to  do  next  is  to  derive  an  expression  for  At'  to  be  substituted  into  (3.3-102),  (3.3- 
108),  and  (3.3-103)  so  that  the  unit  vectors  given  by  (3.3-100)  and  (3.3-101)  can  be  computed. 

The  time  delay  x  (the  amount  of  time  it  takes  for  the  transmitted  acoustic  signal  to  begin  to 
appear  at  the  receiver)  can  be  expressed  as  (see  Fig.  3.3-3) 


x  = 


c  c 


(3.3-109) 


Since 


r'A  =  cAt', 


(3.3-110) 


substituting  (3.3-61),  (3.3-108),  and  (3.3-110)  into  (3.3-109)  yields 


(3.3-111) 


where  x  is  given  by  (3.3-95),  ru2  is  given  by  (3.3-98),  and  the  three  beta  parameters  are  given  by 

(3.3-90)  through  (3.3-92).  Equation  (3.3-1 1 1)  is  the  desired  expression  for  At'  that  will  allow  us  to 
compute  the  unit  vectors  according  to  (3.3-100)  through  (3.3-105)  and  (3.3-108).  Note  that  (3.3- 
111)  can  also  be  obtained  directly  from  (3.3-71)  as  follows.  If  time  t  in  (3.3-71)  is  set  equal  to  the 
time  instant  when  the  scattered  acoustic  field  begins  to  appear  at  the  receiver,  that  is,  if  t  —  x ,  then 
(3.3-71)  reduces  to  (3.3-111).  The  parameter  At'  given  by  (3.3-111)  is  a  constant  whereas  At' 
given  by  (3.3-71)  is  a  function  of  time  t  where  t  >  x . 

The  general  relationship  between  the  acoustic  pressure  p(t,  r)  in  pascals  (Pa)  and  the  velocity 
potential  yM(t,  r)  in  squared-meters  per  second  is  given  by 


(3.3-112) 


where  p0(r)  is  the  ambient  (equilibrium)  density  of  the  ocean  in  kilograms  per  cubic  meter,  shown 

as  a  function  of  position  r  =  (x,y,z).  Note  that  lPa  =  lN/m2  where  1 N  =  1  kg -m /sec2. 
Substituting  (3.3-87)  into  (3.3-86)  yields 
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iOTt  'o,l  M, 2 


t>T. 


(3.3-113) 


Substituting  (3.3-113)  into  (3.3-112),  and  since  the  unit  vectors  n'  ,  and  h"2  are  constants,  and  the 
ambient  density  is  constant  in  our  problem,  the  acoustic  pressure  in  pascals  incident  upon  the 
receiver  at  r2"  =  (x",y",z")  is  given  by 


pit,  r")  =  -j27tfsp0S0HM(t,r''\f,r0)exp(+j2nft),  t>r. 


(3.3-114) 


where  HM (t, r2"  | /,  r0 )  is  given  by  (3.3-87). 

Let  us  next  relate  the  scattering  amplitude  function,  the  differential  scattering  cross  section, 
and  target  strength  of  the  discrete  point  scatterer.  The  target  strength  (TS)  is  defined  as  follows 
[13]: 


TS  =  10  log  10 


/  r  y\/  r,  ff  \ 

A, 


1ref 


dB  re  A 


ref  ’ 


(3.3-115) 


where  [12,  13] 


(3.3-116) 


is  the  differential  scattering  cross  section  with  units  of  squared  meters,  Img.(  r/)  and  Imgi  ( r2  )  are 
the  time-average,  incident  and  scattered  intensities,  respectively,  with  units  of  watts  per  squared 
meter,  gx(f,novh”z)  is  the  scattering  amplitude  function  of  the  discrete  point  scatterer  with  units 
of  meters,  and  An,  is  a  reference  cross-sectional  area  commonly  chosen  to  be  equal  to  1  m2. 


Example  3.3-1  Monostatic  Scattering  Geometry 

We  begin  this  example  by  considering  the  case  when  all  three  platforms  are  in  motion  with 
constant  velocity  vectors  V0,  V,,  and  V2,  but  the  scattering  geometry  is  monostatic  versus  bistatic. 
For  a  monostatic  (backscatter)  scattering  geometry,  both  the  transmitter  and  receiver  are  located  at 
the  same  position,  that  is, 


r2  =  r0. 


Therefore  (see  Fig.  3.3-1), 


v2  =  v0, 


(3.3-117) 


(3.3-118) 
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(3.3-119) 


ri,2=r0.P 


and 


«u  =  -«o.i-  (3.3-120) 

With  the  use  of  (3.3-118)  through  (3.3-120),  the  time-variant,  space-variant,  complex  frequency 
response  of  the  ocean  given  by  (3.3-87)  reduces  to 

HM(t,r"  |/,r0)  *  g1(f,n'0>l,hf2)  1  °'^  exp[+;'2^/5  (t  -  T)]txp{-jljzft), 

(4  7tr0l) 

(3.3-121) 

where  r"  is  given  by  (3.3-99)  and  (3.3-111),  n'0  l  is  given  by  (3.3-100),  (3.3-102),  and  (3.3-111), 
n"  is  given  by  (3.3-101),  (3.3-103),  (3.3-108),  and  (3.3-111),  jSg”  = -jSg0 ,  pff  =  -J8jf , 

P  1,2  Po.l  ♦ 


[i-/C] 

[l  +  j8Sn] 

[i+C”]-| 

is  the  dimensionless,  time-compression  /  time-stretch  factor, 

^(1.0)  _  n0,l  *  X,o 
0,1  c 

is  the  dimensionless,  normalized  radial-component  of  V1>0  in  the  direction  n01, 

0(2,1)  _  «0.1  *  Vw 

Po.i  — 

c 

is  the  dimensionless,  normalized  radial-component  of  V21  in  the  direction  n01, 


(3.3-122) 


(3.3-123) 


(3.3-124) 


0(2)  _  «0,1#V2  _«0,1*V0 
Po.l  —  — 


(3.3-125) 


is  the  dimensionless,  normalized  radial-component  of  V2  =  V0  in  the  direction  n01, 

Vlt0  =  V,  -  (nVl  •  V0)nv,  =  [V;  -  V0(^  •  nVo )]«,,  (3.3-126) 

is  the  velocity  vector  of  the  discrete  point  scatterer  relative  to  the  velocity  vector  of  the  transmitter  in 
the  direction  of  the  velocity  vector  of  the  discrete  point  scatterer  hVi , 
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(3.3-127) 


V2,  =  V2  -  (nVt  •  V,)nVi  =  [v,  -  V,(nv>  )]«* 

-  V0 - («v„ •  V, )«„,  =[%-!' («,, . < )] 

is  the  velocity  vector  of  the  receiver  relative  to  the  velocity  vector  of  the  discrete  point  scatterer  in  the 
direction  of  the  velocity  vector  of  the  receiver  hv^  -  hVo , 


T 


e _ 1  _  /P(2,D  „ 

^  1  Hqa  'o.l 

i+flS”  c 


(3.3-128) 


is  the  time  delay  in  seconds, 


h  = 


2ro,i 


(3.3-129) 


is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmitter  and  receiver 
are  in  motion,  and  r0 ,  is  given  by  (3.3-97). 

If  both  the  transmitter  and  receiver  are  not  in  motion,  then 

V2=V0=0.  (3.3-130) 

Therefore, 

V^Vp  (3-3-131) 

and  since  the  relative  velocity  vector  V21  given  by  (3.3-127)  is  undefined  when  V2  =  V0  =  0,  we  set 


In  addition  (see  Fig.  3.3-3), 


vw  =  o. 


ri,2=ri.2=r0,l 


(3.3-132) 


(3.3-133) 


and 


and  from  (3.3-99), 


A//  A/  /V/ 

^1,2  —  ^1,2  ^ 


r/r  —  r 
l2  h- 


'0,1’ 


(3.3-134) 


(3.3-135) 


With  the  use  of  (3.3-130)  through  (3.3-132),  (3.3-134),  and  (3.3-135),  the  time-variant,  space- 
variant,  complex  frequency  response  of  the  ocean  given  by  (3.3-121)  reduces  to 

HM(t,r2  |/,r0)  =  *,( /,^.p-<i)^^5^ex $+j2nfs(t  -  x)]exp{-j2nft), 

v*™  r0 1 ) 


(3.3-136) 
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where  n'0l  is  given  by  (3.3-100),  (3.3-102),  and  (3.3-111), 

0(1,°)  _  _  Od.O)  =  _  0(1)  0(2)  _  o(2)  =  q 

A-*  1,2  A>0,1  K0.1’  A-' 1,2  Po,l  u’ 


/C=/r 


0,1’ 


(2,1) 


A, 


o. 


i+flK 

(3.3-137) 

is  the  dimensionless,  time-compression  /  time-stretch  factor. 

P 

72  •  V 

'(1)  _  "0,1  M 

0,1  “ 

c 

(3.3-138) 

is  the  dimensionless,  normalized  radial-component  of  V,  in  the  direction  n0 ,, 

1 

x  —  — 
s 

f  1 

l— H 

1 

03 

s — ✓ 

+ 

1 _ 1 

(3.3-139) 

is  the  time  delay  in  seconds. 

_  2ro,i 

To  = 

c 

(3.3-140) 

is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmitter  and  receiver 
are  in  motion,  and  r0 ,  is  given  by  (3.3-97).  Note  that  (3.3-136)  through  (3.3-140)  are  identical  with 

(3.2-105)  through  (3.2-109),  respectively,  derived  in  Example  3.2-1. 

Finally,  if  in  addition  to  both  the  transmitter  and  receiver  not  being  in  motion,  the  discrete 

point  scatterer  is  also  not  in  motion,  then 

v,  =  o. 

(3.3-141) 

Therefore, 

s  =  1, 

(3.3-142) 

and  from  Fig.  3.3-3, 

<i  =  ro,i 

(3.3-143) 

and 

«0,1=«0,I- 

(3.3-144) 

With  the  use  of  (3.3-142)  and  (3.3-144),  the  time-variant ,  space-variant,  complex  frequency 
response  of  the  ocean  given  by  (3.3-136)  reduces  to  the  following  exact  time-invariant,  space- 
variant,  complex  frequency  response  of  the  ocean: 
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/  \  /  \  /  >  exp[-2a(/)r0 ,1 

HM(t,r2  |/,r0)  =  HM(f,r2\r0)  =  8l(f,n 0A,-n0A) - l(4  ^■—■Qxp(-j2nft), 


where  n01  is  given  by  (3.1-18), 


(3.3-145) 


(3.3-146) 


is  the  time  delay  in  seconds,  and  r0A  is  given  by  (3.3-97).  Note  that  (3.3-145)  and  (3.3-146)  are 
identical  with  (3.1-32)  and  (3.1-33),  respectively,  derived  in  Example  3.1-1,  and  with  (3.2-1 18)  and 
(3.2-1 19),  respectively,  derived  in  Example  3.2-2. 


Example  3.3-2  Bistatic  Scattering  Geometry 

We  begin  this  example  by  considering  the  case  when  only  the  discrete  point  scatterer  is  in 
motion  with  constant  velocity  vector  V,  and  the  scattering  geometry  is  bistatic.  Therefore, 


Vo  =  0 


(3.3-147) 


v2  =  o. 


(3.3-148) 


and  as  a  result, 


V  =  V 

vl,0  VP 


(3.3-149) 


and  since  the  relative  velocity  vector  V2 ,  given  by  (3.3-94)  is  undefined  when  V2  —  0,  we  set 


v2,  =  o. 


(3.3-150) 


In  addition  (see  Fig.  3.3-3), 


r"  =  r' 
'1,2  '1,2 


(3.3-151) 


«1,2  =  «1,2’ 


(3.3-152) 


and  from  (3.3-99), 


r2  =r2. 


(3.3-153) 


With  the  use  of  (3.3-148)  through  (3.3-150),  (3.3-152),  and  (3.3-153),  the  time-variant,  space- 
variant,  complex  frequency  response  of  the  ocean  given  by  (3.3-87)  reduces  to 
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H 


(t,r2  |/,r0)  -  gx(f,h o,Pn;,2) £XP^  ^)(r°-1-+  !±S. exp[+j2jifs(t  -  z)]exp(-j2nft), 

Lo7i  r0  [  rl2 


(3.3-154) 

where  n0 ,  is  given  by  (3.2-85),  (3.2-87),  and  (3.2-94),  n,'2  is  given  by  (3.2-86),  (3.2-88),  and  (3.2- 
94), 


s  = 


1-C 


1 

is  the  dimensionless,  time-compression  /  time-stretch  factor, 

pw=Vl L  (3.3-156) 

c 

is  the  dimensionless,  normalized  radial-component  of  V,  in  the  direction  n0 ,, 

(3.3-157) 

c 

is  the  dimensionless,  normalized  radial-component  of  V,  in  the  direction  h12. 


(3.3-155) 


T 


1 

5 


T0  +  (J-1) 


is  the  time  delay  in  seconds. 


*o  = 


r0A  +  rU2 


(3.3-158) 


(3.3-159) 


is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmitter  and  receiver 
are  in  motion,  and  r0 ,  and  r, ,  are  given  by  (3.3-97)  and  (3.3-98),  respectively.  Note  that  (3.3-154) 

through  (3.3-159)  are  identical  with  (3.2-77)  through  (3.2-82),  respectively. 

Finally,  if  the  discrete  point  scatterer  is  not  in  motion,  then 


V,=0. 

(3.3-160) 

Therefore, 

5  =  1, 

(3.3-161) 

and  from  Fig.  3.3-3, 

ro,i  =  ro,v 

(3.3-162) 

«0.1  =  «0.1  * 
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(3.3-163) 

(3.3-164) 


ri,2  =  ri,n 


and 


(3.3-165) 


With  the  use  of  (3.3-161),  (3.3-163),  and  (3.3-165),  the  time-variant,  space-variant,  complex 
frequency  response  of  the  ocean  given  by  (3.3-154)  reduces  to  the  following  exact  time-invariant, 
space- variant,  complex  frequency  response  of  the  ocean: 


/  ,  \  /  |  \  A  \  exp[~ a(f)(  ro,i +  ri,2 )] 

2  |/’ro)  -  ^w(/’r2  lro)  ~  £i(/’no,i’ni,2j  16tr2r01rl2 


txp(-j2nfr). 


(3.3-166) 


where  h0 ,  and  hl  2  are  given  by  (3.1-18)  and  (3.1-19),  respectively, 

T=  r<M  +  yi,2  (3.3-167) 

c 

is  the  time  delay  in  seconds,  and  rw  and  r12  are  given  by  (3.3-97)  and  (3.3-98),  respectively.  Note 
that  (3.3-166)  and  (3.3-167)  are  identical  with  (3.1-17)  and  (3.1-20),  respectively,  and  with  (3.2- 
1 16)  and  (3.2-1 17),  respectively,  derived  in  Example  3.2-2. 


Example  3.3-3  Pulse  Propagation 

In  this  example  we  will  demonstrate  how  to  use  the  pulse-propagation  coupling  equations 
given  by  (2.2-44)  through  (2.2-47)  for  the  bistatic  scattering  problem  shown  in  Fig.  3.3-1.  We 
begin  by  evaluating  (2.2-47),  which  is  repeated  below  for  convenience: 

yw(t,r" )  =  E  JI WK  (/  +  fc)STXf  +  fc)HM(t,r'R[  \f  +  fc,rTi  yxp(+j2nft)df , 

i~l 

(2.2-47) 


where  yM(t, r"  )  is  the  complex  envelope  of  the  acoustic  signal  (velocity  potential)  incident  upon 
element  n  in  the  receive  array  at  time  t  and  position  =  {Xr^Jr^Zr)  with  units  of  squared- 
meters  per  second. 

The  time-variant,  space-variant,  complex  frequency  response  of  the  ocean  that  describes  the 
bistatic  scattering  problem  shown  in  Fig.  3.3-1  is  given  by  (3.3-87).  Therefore 


xexp[-j2n(f  +  fc)t]. 
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(3.3-168) 


where 


[i-ft:?] 

1-  #(2,1) 

1  P  i.«„J 

1- 

_  0(1,0)  1 
Pl.*.\ 

+ 

is  the  dimensionless,  time-compression  /  time-stretch  factor  [see  (3.3-88)], 

o(l,0)  _  nT„\  *  ^1,0 

Pt,  4  “ 


(3.3-169) 


(3.3-170) 


is  the  dimensionless,  normalized  radial-component  of  V,  0  in  the  direction  hT  X  [see  (3.3-89)  and 
(3.1-37)], 


P 


(2,1) 

1,R„ 


n 


1,R 


•V, 


2,1 


C 


(3.3-171) 


is  the  dimensionless,  normalized  radial-component  of  V21  in  the  direction  nlR^  [see  (3.3-90)  and 
(3.1-38)], 


0(1.0)  _  nUR. 

Pi.*.-' 


•vw 


(3.3-172) 


is  the  dimensionless,  normalized  radial-component  of  V,  0  in  the  direction  ^  [see  (3.3-91)  and 
(3.1-38)], 


c 


(3.3-173) 


is  the  dimensionless,  normalized  radial-component  of  V2  in  the  direction  hlR^  [see  (3.3-92)  and 
(3.1-38)], 


v,.„  =  V,  -  (nv<  •  V0)Bv,  =  [15  -  V0(nVi  •  K, )]«, 


(3.3-174) 


is  the  velocity  vector  of  the  discrete  point  scatterer  relative  to  the  velocity  vector  of  the  transmitter  in 
the  direction  of  the  velocity  vector  of  the  discrete  point  scatterer  hVi  [see  (3.3-93)], 


v,,  =  V,  -  (K,  • V,  K,  =  [v2  -  V,  (n*  •  nv,  )]«* 


(3.3-175) 


is  the  velocity  vector  of  the  receiver  relative  to  the  velocity  vector  of  the  discrete  point  scatterer  in  the 
direction  of  the  velocity  vector  of  the  receiver  nVj  [see  (3.3-94)], 


ZTi,R~' 


%*» 


*0  + 


ST:,R„  ^  Pi ,R„  ri.R 

1  -ffil  c 


(3.3-176) 


87 


is  the  time  delay  in  seconds  (the  amount  of  time  it  takes  for  the  acoustic  signal  transmitted  from 
element  i  in  the  transmit  array  to  begin  to  appear  at  element  n  in  the  receive  array)  [see  (3.3-95)], 


T 


0  — 


r7]  .l 

c 


(3.3-177) 


is  the  time  delay  in  seconds  when  neither  the  discrete  point  scatterer  nor  the  transmit  and  receive 
arrays  are  in  motion  [see  (3.3-96)],  and 


rT,,\  ~ 


(3.3-178) 


and 


ri,R„  ~ 


(3.3-179) 


are  the  distances  (ranges)  in  meters  measured  from  element  i  in  the  transmit  array  to  the  discrete 
point  scatterer,  and  from  the  discrete  point  scatterer  to  element  n  in  the  receive  array,  respectively 
[see  (3.3-97)  and  (3.3-98)],  when  motion  begins  at  time  t  =  t0  =  0  seconds. 

The  position  vector  from  the  origin  to  element  n  in  the  receive  array  when  the  scattered 
acoustic  field  is  incident  upon  element  n  is  given  by  [see  (3.3-99)] 

tr„  =  rR„  +  At'V2  +  (t-  At' )  V2>1 ,  t>TT^K,  (3.3-180) 


where  rR  is  the  position  vector  from  the  origin  to  element  n  in  the  receive  array  when  motion 
begins  at  time  t  =  t0  =  0  seconds,  and  [see  (3.3-71)] 


At'~ 


t>T 


T,,R„ 


(3.3-181) 


In  addition,  in  order  to  evaluate  the  scattering  amplitude  function,  we  need  the  following 
equations: 


A/ 

n 


_  rrf.i 


7},1  1 


r„i| 


(3.3-182) 


is  the  dimensionless  unit  vector  in  the  direction  measured  from  element  i  in  the  transmit  array  to  die 
discrete  point  scatterer  when  the  transmitted  acoustic  field  is  first  incident  upon  the  discrete  point 
scatterer  at  time  t '  [see  (3.3-100)],  and 


ll,R„ 


1,R„ 


(3.3-183) 


is  the  dimensionless  unit  vector  in  the  direction  measured  from  the  discrete  point  scatterer  to 
element  n  in  the  receive  array  when  the  scattered  acoustic  field  is  first  incident  upon  element  n  [see 
(3.3-101)],  where  [see  (3.3-102)] 
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rr  a  —  rr  .i +  ^  vi.o 

(3.3-184) 

and  [see  (3.3-103)] 

KR=rUR+An\2-\0)  +  At"Y2A, 

(3.3-185) 

where  [see  (3.1-4)] 

rr„i  =  ri  -  rr; 

(3.3-186) 

and  [see  (3.1-10)] 

r,  - =r, -r, 

(3.3-187) 

are  the  position  vectors  measured  from  element  i  in  the  transmit  array  to  the  discrete  point  scatterer, 
and  from  the  discrete  point  scatterer  to  element  n  in  the  receive  array,  respectively,  when  motion 
begins  at  time  t  =  t0  =  0  seconds,  where  [see  (3.3-1 1 1)] 


A  t'~ 


i  -  ff£l+ ffi.- 


and  [see  (3.3-108)] 


&t"  =  TTjR-  At' , 


(3.3-188) 


(3.3-189) 


where  At'  is  given  by  (3.3-188).  Note  that  the  parameter  At'  given  by  (3.3-181)  is  a  function  of 
time  t  whereas  At'  given  by  (3.3-188)  depends  on  the  time-delay  Tr  ^ . 

Substituting  (3.3-168)  into  (2.2-47)  yields 


Sw )  : 


ej: 

i=  1 


^ ( / )  exp[+ finf- st,.rS*  ~  TT,.*j\dfexp[+j27rfc  si i.K.(f  ~  tt )]exp(-;2^/f  0, 


(3.3-190) 


or 


yAf(t’rR„  )  ~  Xri,R„(Srl.R„[t  ^Ti,R])eXp\+J2ftfc(STi,R„  l)t]eXp(  j^fc  ST„Rm^T„R.^> 

(3.3-191) 


i=i 


where 


/  Nexp[  -a(f  +  fc)(rTl  +  rlR) 

Z'^(/)  =  X(/)c,(/  +  /c)S,(/  +  /c)gl  (/  +  /e,  - lA 

Tf,\  1  ,R„ 


(3.3-192) 
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and 


4„»,«)  =  F/'{^.,.(/)}  =  £^,„,(/)exp(+y2^/0#.  (3.3-193) 

Equation  (3.3-191)  indicates  that  the  complex  envelope,  yM(t, r" ),  of  the  acoustic  signal  (velocity 

potential)  incident  upon  element  n  in  the  receive  array  is  equal  to  the  sum  of  time-compressed  / 
time-stretched,  time-delayed,  Doppler  shifted,  modified  versions  of  the  complex  envelopes  of  the 
transmitted  electrical  signals  from  each  element  in  the  transmit  array.  The  amount  of  Doppler  shift 
in  hertz  is  given  by 

*r, <33-194> 

Equation  (3.3-192)  shows  how  the  frequency  spectrum  of  the  complex  envelope  of  the  transmitted 
electrical  signal,  X(f),  gets  modified  by  various  other  important  frequency-dependent  functions  in 
order  to  produce  the  complex  envelope,  yM(t, r*' ),  of  the  acoustic  signal  (velocity  potential)  incident 

upon  element  n  in  the  receive  array. 

Before  continuing  with  the  analysis,  several  comments  concerning  the  dimensionless,  time- 
compression  /  time-stretch  factor  sTRn  and  (3.3-194)  are  in  order.  There  are  three  ranges  of  values 

for  sT  R  ;  namely,  sT.  <R=  1,  sTi  R>  1,  and  sTR<  1.  When  sT,R  =  1,  this  is  an  indication  that  there  is 

no  time  compression  or  time  stretch,  and  as  a  result,  there  is  no  change  in  signal  bandwidth.  Also, 
when  sT  R  =1,  <t>TR  =  0;  that  is,  there  is  no  Doppler  shift  [see  (3.3-194)].  For  example,  if  there  is 

no  motion,  then  all  three  velocity  vectors  are  equal  to  the  zero  vector.  Therefore,  all  four  beta 
parameters  are  equal  to  zero  and,  as  a  result,  sT  R  =  1 .  When  sTi  R^>  1,  this  is  an  indication  of  time 

compression,  that  is,  the  duration  of  the  received  signal  given  by  (3.3-191)  is  decreased  (compared 
to  the  duration  of  the  transmitted  signal)  resulting  in  an  increase  in  signal  bandwidth.  Also,  when 
sT  R  >1,  <pT  R  >  0 ;  which  indicates  a  positive  Doppler  shift  [see  (3.3-194)].  Finally,  when  sT,  R<  1, 

this  is  an  indication  of  time  stretch,  that  is,  the  duration  of  the  received  signal  is  increased 
(compared  to  the  duration  of  the  transmitted  signal)  resulting  in  a  decrease  in  signal  bandwidth. 
Also,  when  sT  R  <1,  <pTR  <  0 ;  which  indicates  a  negative  Doppler  shift  [see  (3.3-194)]. 

If  we  next  substitute  (3.3-191)  into  (2.2-46),  then  we  obtain 

)  =  5  - ^.^])exp[+y2^/c sT^p-  tT^ )]}.  (3.3-195) 

1=1  L 

Since  x(t)  given  by  (2.2-28)  and  x(t)  given  by  (2.2-30)  are  related  by 

x(t)  =  Re{ic(t)exp(+ j2nfct)},  (3.3-196) 

if  we  let 

x'Tj<R(t)  =  Re jx'ri  R{t) txp(+j27tfc  r)} ,  (3.3-197) 

then  (3.3-195)  can  be  rewritten  as 
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(3.3-198) 


T 

y^rRn )  ~  ^  —  ^V|f /?„])* 


i=l 


Substituting  (3.3-198)  into  (2.2-45)  yields 


1 


)  ~  |  i^7;^fl(7?/*y7;.,^)exp(  j2nr]TTiR),  (3.3-199) 


i=l  15. 


T,/?r 


and  upon  substituting  (3.3-199)  into  (2.2-44),  we  obtain 


y(t,r")  *  2  £  X",K  mxV[+jlnr)(t  -  rTi,R)\dr),  n  =  l,2,...,iV*,  (3.3-200) 

1=1 


or,  finally. 


y{^zR)  ~  ^  xT_  R(t  —  zTR),  n  =  1,2  ,...,Nr, 


(3.3-201) 


1=1 


where 


and 


^Ti.R^V)  ~  I  I  ^r;  ,R„(77/ ST>,R„  )CR„  (v) Sr, (y) 


*Ti,Rn 


(3.3-202) 


*W*>  =  F;'{X".  „,(!))}  =  £  X'',„,(i;)exp(+y2®l)()<il). 


(3.3-203) 


Equation  (3.3-201)  indicates  that  the  output  electrical  signal  y(t, rR)  from  element  n  in  the  receive 

array  at  time  t  and  position  rR  =  (xR  ,yRn,zRJ  with  units  of  volts  is  equal  to  the  sum  of  time- 

delayed,  modified  versions  of  the  transmitted  electrical  signals  from  each  element  in  the  transmit 
array. 
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4  Summary 

A  set  of  pulse-propagation  coupling  equations  was  successfully  derived.  They  couple  the 
output  electrical  signal  at  a  point  element  in  a  receive  array  to  the  transmitted  electrical  signal  at  the 
input  to  a  transmit  array  via  the  complex  frequency  response  of  a  fluid  medium  (e.g.,  air  or  water). 
The  pulse-propagation  coupling  equations  are  based  on  linear,  time-variant,  space-variant,  filter 
theory,  the  principles  of  complex  aperture  theory  and  array  theory,  and  solving  a  linear  wave 
equation,  which  includes  satisfying  all  boundary  conditions,  including  the  boundary  condition  at  the 
source.  The  time- variant,  space- variant,  complex  frequency  response  of  the  ocean  was  shown  to  be 
the  time-harmonic  solution  of  a  linear  wave  equation  when  the  source  distribution  is  a  time- 
harmonic  point  source. 

The  pulse-propagation  coupling  equations  provide  a  consistent,  logical,  and  straightforward 
mathematical  framework  that  can  be  used  to  accurately  model  the  propagation  of  small-amplitude 
acoustic  pulses  in  the  ocean  for  a  bistatic  scattering  problem.  The  main  features  of  the  pulse- 
propagation  coupling  equations  are  as  follows:  1)  transmitted  electrical  signals  are  modeled  as 
amplitude-and-angle-modulated  carriers,  2)  both  the  transmit  and  receive  apertures  are  modeled  as 
volume,  conformal  arrays  composed  of  unevenly-spaced,  complex-weighted,  point  elements  (this 
type  of  model  for  both  of  the  apertures  allows  for  maximum  flexibility),  3)  the  complex  weights  are 
frequency  dependent  and  allow  for  beamforming,  4)  the  performance  of  the  point  elements  in  both 
the  transmit  and  receive  arrays  are  characterized  by  frequency-dependent,  transmitter  and  receiver 
sensitivity  functions,  and  5)  the  solution  of  a  linear  wave  equation  is  given  by  the  complex 
frequency  response  of  the  fluid  medium.  It  is  important  to  note  that  attention  to  all  proper  units  of 
measurement  were  taken  into  account  in  order  to  ensure  the  accurate  prediction  of  signal  strength 
levels  at  each  element  in  a  receive  array.  This  is  especially  important,  for  example,  in  order  to  obtain 
accurate  probability  of  detection  results. 

Derivations  of  the  complex  frequency  response  of  the  ocean  for  the  following  three  different 
bistatic  scattering  problems  were  performed:  1)  no  motion,  2)  only  the  discrete  point  scatterer  is  in 
motion,  and  3)  all  three  platforms  (the  transmitter,  discrete  point  scatterer,  and  receiver)  are  in 
motion.  Specific  examples  on  the  use  of  the  pulse-propagation  coupling  equations  were  given  for 
the  three  different  bistatic  scattering  problems.  Scatter  from  a  discrete  point  scatterer  was  modeled 
via  the  scattering  amplitude  function,  which  is  a  complex  function  (magnitude  and  phase)  and  is,  in 
general,  a  function  of  frequency,  the  direction  of  wave  propagation  from  the  source  to  the  scatterer, 
and  the  direction  of  wave  propagation  from  the  scatterer  to  the  receiver.  In  addition  to  the  scattering 
amplitude  function,  frequency-dependent  attenuation  was  taken  into  account  in  order  to  model  the 
propagation  of  sound  from  transmitter  to  discrete  point  scatterer,  and  from  discrete  point  scatterer  to 
receiver.  The  dimensionless,  time-compression  /  time-stretch  factor  was  derived  and  discussed  for 
the  two  bistatic  scattering  problems  involving  motion.  The  time-compression  /  time-stretch  factor 
takes  into  account  the  relativistic  effects  of  motion  and  provides  for  accurate  time  delay  and  Doppler 
shift  values. 
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